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REDUCE II - A COMPUTER PROGRAM FOR THE SYMBOLIC REDUCTION OF 

LARGE BLOCK DIAGRAMS 
by Carl F. Lorenzo and John P. Riehl 
Lewis Research Center 

SUMMARY 

A computer program, REDUCE II, has been written to reduce very large (linear) 
block diagrams. The program requires as input a set of algebraic expressions repre- 
senting the block diagram, the desired transfer function(s), and a string of variables 
indicating the desired order of reduction. REDUCE II then calculates the transfer func- 
tions) of any output variable to any or all input variables. The program is premised on 
reduction, by loops, of the block diagram. The solution is presented in the compact form 
of a set of nested functions (super G’s), The super G f s define the transfer functions in 
terms of the system information (G T s) and the preceding super G’s, 

The program can handle as many as 600 equations and 600 variables and is intended 
as a tool for the analysis of complex control and dynamic systems. Several applications 
are presented. 

A FORTRAN evaluation program, EVAL II, accepts the symbolic output generated 
by REDUCE D and evaluates it numerically to obtain amplitude ratio and phase angle as 
functions of frequency for any desired G functions. 


INTRODUCTION 

The use of block diagrams for the analysis of dynamic problems is well known and 
widespread. As the systems of interest increase in complexity, new techniques are 
sought to aid the controls engineer in the analysis and reduction of such systems. 

Indeed, at least part of the motivation of modern control theory is to deal with the 
problem of large system complexities through the organizational features provided by 
the matrix approach. 

The concept of block diagram reduction is to simplify the complex diagram to an 
equivalent diagram composed of a single block or a few blocks, that is, to topologically 
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simplify the block diagram. 

In reference 1, the predecessor of this effort, a program was evolved which sym- 
bolically reduced block diagrams to obtain S form and io> form algebraic solutions for 
desired transfer functions. This approach suffices for intermediate- size diagrams and 
also gives fully expanded closed-form solutions. However, the method is limited to sys- 
tems of approximately 30 equations (for transfer functions (G’s) which are not too com- 
plicated) due to storage limitations (for the IBM 7094 computer). Experience with the 
predecessor program has shown that while the closed-form solutions are desirable for 
many scientific endeavors, the bulk of the potential applications were of an engineering 
nature and could not be performed due to the size limitation. Most of these engineering 
applications required only numerical results representing the desired transfer function. 
However, the applications tended to be large and quite complex. 

Direct matrix methods (ref. 2), for example, could be applied to these systems. 
However, these methods would require a complete matrix manipulation for each frequen- 
cy and would become quite expensive in terms of computer time. 

Another possible approach is to use the symbolic techniques as a part of a reduction 
scheme, together with a decompositional method which would not allow the solution to 
grow within the computer. The development of this approach is the objective of this 
study. 

The report develops the mathematical concepts and a computer program (REDUCE 
II) to implement this approach. A numerical evaluation program (EVAL II) to use these 
symbolic results is also developed. Applications of the programs are presented. 


REDUCTION TECHNIQUE 

In general, dynamic and control systems are analyzed by considering the character- 
istics of the elements (individual transfer functions) and combining these elements into a 
system which is the block diagram. The block diagram elements are (1) summers, 
which add or subtract two or more signals X’S; (2) blocks, which accept one signal and 
modify it through multiplication by a frequency-sensitive operator G(S); and (3) nodes, 
which split a signal into two or more parts. These elements are then combined to yield 
linear block diagrams. For details of block diagrams and their reductions, refer to ref- 
erence 1. This study addresses itself to the reduction of very large linear block dia- 
grams (of arbitrary topology) to obtain the transfer functions of a desired output to any 
or all inputs. The symbolic computer language FORMAC (ref. 3) is used to achieve the 
reduction. The FORMAC language is an experimental extension of FORTRAN IV for the 
IBM 7094 computer. 
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Basic Philosophy 


The fundamental philosophy of reduction used in the REDUCE II program entails two 
points which are incorporated to maximize block-diagram-size capability. First, any 
steps which promote significant expansion of the forms within the computer are avoided. 
For example, substitution of values of the G's is not done in the reduction program. 
Second, the reduction of the block diagram must proceed in such a manner as to have the 
system of equations collapse rather than grow in the computer. That is, as variables 
are eliminated in the reduction process, the coefficients of the remaining variables be- 
come combinations of the original coefficients (G's). As this occurs, these combinations 
are replaced by single- element coefficients (super G's), and the mathematical equation 
defined by this substitution is immediately output from the computer. 


Technique 

The general approach to the reduction is as follows : 

(1) The system of linear algebraic equations which represent the block diagram is 
input to the computer in symbolic form. 

(2) Input is an order string, that is, a string of variables which guides the manner 
in which the reduction is to be performed and indicates when a substitution is to be 
made. From these data, the program reduces the block diagram by loops and, after 
each loop is reduced, substitutes a new G (a super G) for combinations of G’s formed 
by the reduction process. By reduction of a loop is meant the elimination of those vari- 
ables which form a cycle or a circuit in the graph theoretical sense (e.g. , see refs. 4 
and 5). By virtue of this technique, the system of equations in the computer can only 
collapse; and hence, the largest system which can be handled is approximately the larg- 
est system of equations which can be initially input to the computer, or in practical 
terms, that set of equations which can be held in core. 

From the mathematical point of view then, we are starting with a set of linear equa- 
tions in the variable X with coefficients (G’s) which are considered to be constant. 
Through the reduction and substitution process, we are transforming the original set of 
equations into a new set of equations which define the super G's in terms of the original 
G’s and previous super G's and into a single equation involving the X's of the transfer 
functions. 

Even within this framework, a large number of ways exist to proceed, since the 
form of the equations composing the solution set has not been stipulated. However, the 
following desirable properties of the solution set can be projected from the intended use: 

(1) Since eventually the solution set is to be evaluated numerically for amplitude 
ratio, phase angle, and so forth, the super G’s formed should be nested. That is, a 
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given super G should depend only on the G T s which precede it and on the original G's of 
the block diagram . 

(2) A zero divide occurs when the denominator of a super G is evaluated to be iden- 
tically zero. The solution set should be of such a form that the likelihood of a zero di- 
vide occurring when the solution set is evaluated is minimal. 

(3) The solution set should be of such a form that subsequent equation scaling is not 
necessary. For example, the process could be carried out in such a manner as to allow 
only super G T s of a sum-product form, that is, 

G SUPER = G a G b G c + G d G e G f + • • * etc - 

Symbols are defined in appendix A. Now, while this form would be highly desirable in 
terms of property 2, it leads to a scaling problem. Experience has shown that for cer- 
tain common block diagram topologies (notably some of the ladder networks) the last 
super G, even for relatively small diagrams, can be a function of u> to a very large 
power. This, of course, would only allow evaluation to a very small value of o> without 
some form of scaling. If scaling is to be done, it should be integral with the formation 
of the solution set. 

(4) Finally, it is desired that the solution set be as compact or concise as possible. 
This is a purely practical requirement, since for very large block diagrams the volume 
of punched cards and printout could become unwieldy. 

Properties 2 and 3 can be conflicting requirements. Consideration of properties 2 
and 3 suggests that a most desirable form for the super G's composing the solution set 
would be 


( 1 ) 


or as close to this form as possible. This form is characterized by the fact that, for 
G's of the same order, the super G would be of net order 1 and since the denominator 
has a unity additive term, the possibility of a zero-divide condition (by allowing any G 
to equal zero) is null. The G = 0 special case is considered important since it is the 
fundamental method of simplifying block diagram topology. It appears that super G's 

of the form of equation (1) can always be made to occur if the block diagram is reduced 
by loops. 

This is demonstrated by the following simple example. Figure 1 is a block diagram 
of a loop which is being considered separately from a larger block diagram which con- 
tains it. The defining equations for the embedded loop are 
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X 1 - X 2 - ^ = 0 


Xj - x 4 - x 6 = 0 

x 3 ' °l x 5 = 0 


Xj - g 2 x 4 = 0 


Eliminating the loop by eliminating the variables X^ and Xg yields 

G,G 9 g 9 

x 9 = LJ_ x x 6 

1 + GjGg 1 + G^Gg 


and 


G-, GiG 9 

Xo = - X t + Xg 

3 1 + G jG 2 1 * GjG 2 


where X 2 and Xg are the loop outputs defined in terms of Xj and X g , the loop inputs. 
The coefficients of X 1 and X g are the super G’s associated with the reduction of the 
loop. These super G f s are seen to be of the same form as equation (1). 




Figure 1. - Block diagram of embedded loop. 


The fact that this kind of a solution form is possible, even for a loop embedded in a 
large block diagram, has been used as a fundamental premise in the structure of the 
computer program presented in this report. Not every equation affected in the reduction 
of a loop will contain the 1 + GG + . . . coefficient form. That is, those equations out- 
,ide the loop, but containing loop variables, will not at this stage of the reduction contain 
.ne 1 + GG + . . . form. They might at a later stage. However, postponement of the 
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super G substitution is not advised, since the coefficients can grow unreasonably large, 
defeating the program purpose. Contingencies arise when the loops are very large or 
difficult to identify. These details are considered later. 


Program Flow 

Before considering the flow of the reduction program it will be helpful to consider 
the relation between the reduction and evaluation pregrams. This relation is indicated 
in figure 2. The flow chart is composed of two basic parts: the algebraic manipulations 


Algebraic manipulations Numerical manipulations 



Figure 2. - Relation of algebraic and numerical programs. 


are done in FORMAC, and the numerical calculations are done in the companion FOR- 
TRAN evaluation program. The expansive step of substituting the values for the G 
functions into the solution set is postponed until the evaluation program. Thus, larger 
block diagrams can be reduced. Also flexibility is added since the G’s can be changed 
at will when the numerical evaluations are performed. Hence, it is not necessary to re- 
peat the reduction of a block diagram in order to change the system. The only require- 
ment is that the topology of the block diagram is not changed. Specifically, the topology 
can be simplified by letting G equal zero, but not made more complex. 

The logic and flow of the reduction program are illustrated in the flow chart of fig- 
ure 3. An understanding of most of the program function is best obtained by referring to 
a sample reduction as it would be performed by the program. The system to be reduced 
is shown in block diagram form in figure 4. 
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/Read 
U) Solution form 
(2) Block diagram inputs 
13) Transfer functions desired 

(4) Equation of block diagram 

(5) Order string 



— j — — 

i/Vrite 
' error 


Figure 3. - Reduction program flow chart. 
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Figure 4. - Block diagram of system to be reduced. 


The defining equations for this system are 


X 1 - X 7 ' X 8 = 0 

(2) 

!" X 3" X 6" X 9 =(I 

(3) 

X 4 - GjXj . 0 

(4) 

x 6 - g 2 x 4 = C 

(5) 

Xg-GjXg-O 

(6) 

Xj - g 4 x 4 = 0 

(7) 

X, - g 5 x 5 = 0 

(8) 


The inputs are X l and X 2 , and the desired output is Xg; that is, form the transfer 
functions Xg/Xj and 

first eliminates the inner loop, that is, the variables Xg and Xg. After these are elim- 
inated, a super G substitution is made and the next loop (X 4 , X ? , Xg, and Xg) is elimi- 
nated. The process is continued until all variables but the output and inputs are elimi- 
nated. A reasonable order string for this problem then is 

Xg,Xg, , 

where the single commas denote a variable elimination and the double commas denote a 
variable elimination followed by a super G substitution for thoBe coefficients containing 
more than two G's. A discussion of the order string selection is included in appendix B- 
The system of equations is reduced in the following steps in the same manner as 


X4, Xrj f Xg, Xg, , 


X g /X2. Since the program reduces by loops, the order string 
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would occur in the program based on the order string given: 

STEP 1: Eliminate variable X g . Since Xg occurs in equations (3) and (4), the sim- 
pler expression (eq. (4)) is substituted into the more complex (eq. (3)). Then, multiply- 
ing equation (3) by G 4 and substituting X 4 for G^Xg yields 


- X 4 - GjX 6 - GjX, - 0 

(3a) 

0 = 0 

(4a) 


STEP 2: Eliminate variable X g and make super G substitution. Here X g occurs 
in equations (3a) and (5). Again, results from equation (5) are substituted into equation 
(3a), giving 


Gl X2 - X 4 - G 1 G 2 X 4 - GjX, = 0 

(3b) 

0 = 0 

(5b) 

Combining coefficients yields 


GjXg - (1 + GjG 2 )X 4 - G 1 X 9 = 0 

(3c) 

In making the super G substitution, the coefficients of the X’s 

of all effected equations 

(here only eq. (3c)) are now searched for a ±1 ± < £ TT G > form. If the form occurs, 

the equation is divided through by it and super G’s are formed. 

Hence, equation (3c) 

becomes 


GjJCj ♦ X 4 - G 6 Xj = 0 

(3d) 

where 


- G 1 

G fi = 1 

6 i + QiG 2 



is the first super G. 

STEP 3: Eliminate variable X 4 . Equation (7) is used to eliminate X 4 in equation 
(3d); then 

G 4 G 6 X 2 + X 5 - °4 0 e*9 = 0 (3e) 


0 = 0 


(7a) 
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STEP 4: 

Eliminate variable X^. Here X,y occurs in equations (2) and (8); 

hence, 


X 1 - °5 X 5 - x 8 * 0 

(2a) 


o 

II 

o 

(8a) 

STEP 5: 

Eliminate variable Xg. Using equation (6) in equation (2a) gives 



G 3 X 1 - G 3 G 5 X 5 - Xq ' 0 

(2b) 


0 = 0 

(6a) 


STEP 6: Eliminate variable Xg and make super G substitution. Substituting from 
equation (2b) into equation (3e) gives 

G 4°6 x 2 ♦ *5 - G 3 G 4 G 6 X J + G 3 G 4 G 5 G 6 X 5 = 0 OT 

Combining coefficients yields 

G 3 G 4 G 6 X 1 - G 4 G 6 X 2 - < l ♦ G 3 G 4 G 5 G 6> X 5 = 0 «g) 

Recognizing the 1± form and forming the super G’s gives 

+ GgXg + Xg = 0 (3h) 


where 

G G 3 G 4 G 6 

7 " -d + g 3 g 4 g 5 g 5 ) 

and 

G ~ G 4 G 6 

8 -(1 + G 3 G 4 G 5 Gg) 

This completes the reduction as indicated by the order string. The only remaining 
step is to form the desired transfer functions Xg/Xj and X 5 /X 2 - To form Xg/Xp let 
X 2 = 0 in equation (3h); hence, 

G 7 X 1 * - x 5 
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Likewise for X^/Xg, let X^ = 0 to give 



C 8 X 2 - - X 5 


X5 G 8 



While it did not occur in this reduction, if at some step a variable to be eliminated 
occurs in three or more equations, the simplest equation is used to eliminate it in all 
equations in which it occurs. Also, if a substitution is called for by the order string and 

an effected equation does not have a | ±1 ± £ TT Gj form coefficient, a 1.0 can be add- 

V j k / 

ed to the largest coefficient and this factor used to scale the equation. If this is done, 
the solution is said to be of the ARTIFICIAL form. If 1.0 is not added but the coefficient 
itself is used to scale the equation, the solution is called the NATURAL form. The 
choice of form is made by the user. The details of program implementation are covered 
in appendix C. The descriptions, listings, and flow charts for the computer programs 
are detailed in appendix D. A user’s manual is provided in appendix E. 

Discussion of the companion program EVAL II is also contained in these appendixes. 


APPLICATIONS 

Applications of the program were chosen to achieve two objectives: first, to study 
certain operating characteristics of the program, namely, operating time and size capa- 
bilities; and second, to demonstrate the program for typical controls block diagrams. 


Application 1: Simple Ladder 

For this first study, a simple ladder block diagram of the form of figure 5 was re- 
duced. This topology is typical of finite difference approximations to the wave or diffu- 
sion equations. For wave equations, this would be a string of coupled oscillators. 

The basic equation form for each of the ladder’s 600 rungs is 

Xn + 1 - °A * G n x n + 2 * 0 tor n = 1, 2, . . . , 598 
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Figure 5. - Block diagram of simple ladder network. 


To terminate the ladder, the equation 

X 600 " G 599 X 599 = 0 

was used. 

The reduction order string for this problem is 

X 600’ 5 X 599’ * X 598 5 > ‘ * • X 5’ ,X 4 ,,X 3 

The transfer function desired is that for Xg/X^. Typical input equations and the first 
and last parts of the solution set of equations as generated by the program are presented 
in figure 6. The solution set (fig. 6) is output in the form of FORTRAN expression and 
not as a mathematical equation. That is, mathematically 

G(601) = 

DENOM(l) -G 5Q8 G 599 - 1.0 

The time history of reduction of the 600-equation set is shown in figure 7. As might 
be expected, the slope of the curve indicates that the reduction proceeds more quickly as 
variables are eliminated. The total reduction time for the 600- rung ladder is about 
3500 seconds. 

The repetitive nature of block diagram topology allows the solution to be readily 
checked. Also, the repetition, together with the order of solution chosen, allows one to 
consider the problem as a new ladder one rung shorter after each variable is eliminated 
Hence, it could appear as though 600 ladder networks of lengths 1 to 600 were reduced. 
The initial rate of reduction of each of these ladders against the ladder size is shown in 
figure 8. Here the rates vary from 0. 1 variable per second for the 600-equation ladder 
to 0. 5 variable per second for a 25-equation ladder. 

The 600- equation ladder is the practical upper limit for the computing system used 
for this study. Larger problems could be input; however, peripheral storage (noncore) 
would be required, thus heavily penalizing computing time. 
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THE BLOCK DIAGRAM INPUTS 


* * * OUTPUT FROM REDUCE II * * * 


X(l)$ 

THE SOLUTION IS FOR X(2)$ TO X(l)$ 
THE BLOCK DIAGRAM EXPRESSIONS 


1 

X( 

2) 

G( 

1)*X( 

1)+G( 

l)*x( 

3) $ 

2 

X< 

3) 

G( 

2{*X( 

2)+G( 

2)*X( 

4)$ 

3 

X( 

4) 

G( 

3) *X( 

3)+G{ 

3)*X( 

5)1 

4 

X( 

5) 

G( 

4)*X( 

4)+G( 

4)*X( 

6)$ 

5 

x( 

6) 

G( 

5) *X( 

5)+G( 

5)»X( 

7)* 

6 

X( 

7) 

G( 

6) *X( 

6)+G( 

6)*X( 

8)$ 

7 

X( 

8) 

G( 

7 ) *X( 

7)+G< 

7)*X( 

9)$ 

8 

X( 

9) 

G( 

8)»X< 

8)+G{ 

8)*X( 

10)$ 

9 

X( 

10) 

G( 

9) *x( 

9)+G( 

9)*x( 

11)$ 

10 

X( 

11) 

G( 

10)*X( 

10)+G( 

10)*X( 

12)$ 

11 

x( 

12) 

G( 

ll)*x( 

1 1 ) +G { 

ll)»x( 

13)$ 

12 

X( 

13) 

G( 

12 ) *X( 

12)+G( 

12 ) *X( 

14) $ 

13 

X( 

14) 

G( 

13)*X{ 

13)+G( 

13)*X( 

15)$ 

14 

X( 

15) 

G( 

14)*X( 

14)+g( 

l4)*X( 

16 ) $ 

15 

X( 

16) 

G( 

15) *X( 

15)+G( 

15)*X( 

17)$ 

16 

X( 

17) 

G( 

16) *X( 

l6)+0( 

16)*X( 

18 )$ 

17 

X( 

18) 

G( 

17)*X( 

17)+GC 

17)*X( 

19)$ 

18 

X( 

19) 

G( 

1 8 ) *X t 

18)+G( 

l8)*X( 

20)$ 


593 

X( 

594)-GC 

593)*X( 

593)+G( 

593)*X( 

594 

x( 

595)-G( 

594 ( *X( 

594)+G( 

594 ) *X( 

595 

x( 

596)-G( 

595)*X( 

595)+0( 

595)*X( 

596 

x( 

597)-G( 

596)*X( 

596)+G( 

596) *X( 

597 

x( 

598)-G( 

597) *X< 

597 )+G( 

597)*X( 

598 

x( 

599)-G( 

598) *X( 

598)+G( 

598) *X( 

599 

X( 

600)-G( 

599) *X( 

599)$$ 



DENOM ( 1)~ 

-G(598)*G(599)-I.O 
Q( 600)-G(598) 

G( 600 )=G ( 6 0 0 ) /DENOM ( I) 

DEN0M( 2 )- 

G(59T)*g(600)-1.0 
G( 60I)=G(597) 

G( 6oi)=a( 601)/DENOM( 2) 

DENOM ( 3)- 

G(596)*G(601)-1>0 


DENOM ( 597)** 

G(2)*G(1195)-l-0 
G( 1196)=G(2) 

G(1196)-G( 1196) /DENOM ( 597) 

DENOM ( 598) = 

G(l)*G(1196>-1.0 

595) $ G ( 1197 )=G(l ) 

596) $ 

597) $ G( 1197 )*G( 1197 ) /DENOM ( 598) 

598) $ 

599) $ TRANSFER FUNCTION X( 2)/X( l) 

600 ) $ 

THE NUMERATOR 
0(1197) 

THE DENOMINATOR 
(- 1 * 0 ) 


la) Input expressions. (b) Solution-set equations. 

Figure 6. - Simple-ladder input and solution set 
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Application 2: Complex Ladder 


For this application, a ladder of coupled four-terminal networks was reduced. The 
ladder is composed of 150 rungs, each having two equations; the block diagram is shown 
in figure 9. For a typical application, each rung of such a ladder would represent a dis- 
tributed parameter representation of the wave equation. That is, the transfer functions 
of the blocks would be complex hyperbolic functions or delay times. The topology has 
wide application. 



' v > V V ' ' v 

Rung n = 1 Rung n = 2 Rung n = 150 

Figure 9, - Block diagram of complex ladder network (four-terminal elements). 


The basic equations for each rung of the ladder, for n = 1, 2, . . . , 150, are 

X 2n+1 ‘ G 4n-3 X 2n-l " G 4n-2 X 2n = 0 

and 

X 2n “ °4n-l X 2n+l " G 4n X 2n+2 = 0 

The desired transfer functions are X 2 /X x and X 2 /X 302 . The reduction order is 

X 301’ ' X 300’ X 299’ ’ X 297’ ’ ’ ‘ • » » X 6’ X 5 > * X 4’ X 3* 

The first and last parts of the solution-set input equations are presented in figure 10 

The history of the reduction is shown in figure 11. The total time for the reduction is 

1010 seconds. 

Somewhat larger ladders of this type could be reduced with the program if the di- 
mension of G is increased. This was not attempted in this study. 
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Figure 1L - Reduction history of 300-equation complex 
ladder. 


Application 3: Jet-Engine-Inlet Shock-Position Control 

The block diagram of figure 12 represents a system for control of the position of a 
shock wave in a supersonic inlet for a jet engine using engine fuel flow. "Modern” con- 
trol techniques were used to design a controller to minimize the occurrence of inlet un- 
starts in the presence of a stochastic disturbance. The control designer wished to ob- 
tain specific transfer functions from which system frequency domain behavior could be 
studied. Since a PADE’ approximation to the inlet plant dead time had been made for 
control design purposes, it was important to determine the effect of this approximation 
on dynamic performance. Transfer functions X^/Xp X^/Xg, and X^/Xg were obtained 
for this reason. The equations representing the plant and controller (the block diagram) 
as input to the program are shown in figure 13(a). Also shown is the order of reduction. 

The block diagram inputs are Xp Xg, and Xg. The results of this reduction (us ing 
the ARTIFICIAL form of solution) are given in figure 13(b). From these results, evalu- 
ations with different plants and controllers can be made as long as the topology of the 
block diagram does not become more complex. Essentially, any simplification of the 
topology is possible. 

For example, we can check the plant transfer function by simply evaluating X^/Xg 
with to Ggj equal to zero. This check capability, by specifying certain G's to be 
zero, is useful when working with very large systems. Also, it is nearly always possi- 
ble when the solution is of the ARTIFICIAL form. 

The time required to execute this reduction is 25 seconds. The system is approxi- 
mately equivalent to 50 elementary equations as determined by consideration of the num- 
ber of simple summers and blocks of the block diagram. 
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Figure 12. - Block diagram of jet engine inlet control. 













THE BLOCK DIAGRAM INPUTS 


XI 1)S 
X ( 2 ) S 
X < 3 ) $ 

THE SOLUTION IS FOR X{7>$ TO THE INPUTS 
THE BLOCK OIAGRAM EXPRESSIONS 

1 -XI13) + X(9)-GI 10 >*X I 16)-G( 11 ) »X 1 1 5 ) -Gt 12 > *X( 1* )-G( 131 *X 1 13 ) -Gl 14) *X ( 1 2) -G 11 5 1 »XI 
1 111$ 

2 X(lt+X(6)-X(7)» 

3 X\ 31 -XI 17 ) -X ( 4 ) $ 

4 XI2)*XI4)-Xl5>$ 

5 X ( 6) -G( l ) *X ( 5 } $ 

6 XI7I-XI 18)-X(8)$ 
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THE REDUCT I DM ORDER IS 

1 XII3) tX(li) ,,XI121 »,XI 13) , .XU4I..XI 15). ,X(16) ,.Xt9).X(lB| .X(8)..X(M .X(17),X14) 
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(b) Solution-set equations. 

Figure 13. - Input and solution set for jet-engine-control problem. 
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Figure 13. - Continued. 
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Figure 13. - Concluded. 


A numerical evaluation for the following set of G’s (system information) was made 
by using the EVAL program: 

G = 846. 7(-0. 833x10" 2 S 3 + 0. 1 S 2 - 0. 5 S + 1.0) 

S 6 + 18. 28 S 5 + 155.4 S 4 + 743.9 S 3 + 2038S 2 + 2823S + 846. 7 

G 2 * 846.7 
A S 

G 4 ,G 5 G g = 2.041662, -0.3595, -2.54, 8.723, -1.5025, -146.932 

G 10’ G 11»* * ♦» g 1 5 = 846.7, 2823, 2038, 743.9, 155.4, 18.28 
G 16 ,G 17’ * * -. G 21 = 998,9> 600 - 97 > 155.18, 20.745, 1.300, 0.038 
For the nonapproximate case, the plant transfer function becomes 

n _ _ 7.06 e' 1 * 08 

Gj 

S 3 + 6.33 S 2 + 20.01 S + 7.06 


The frequency range of interest is 0.01 to 4.0 hertz, in steps of 0.01 hertz. (Be- 
cause of time and frequency scaling requirements of the user, frequency was multiplied 
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by 100. ) The execution time required to evaluate these 400 points plus an additional 
600 points for the automated plot routine for all three transfer functions is 32 seconds. 
Plots of typical results for the X 7 /Xj transfer function are presented in figure 14. The 
difference in the plots is of course the effect of the PADE T approximation at the higher 
frequencies . 



lb) Phase angle. 

Figure 14. - Typical transfer function results, for jet-engine-control problem. 


Applications: Some Comments 

Execution time . - A summary plot of reduction rate as a function of problem size 
for the applications studied is presented in figure 15. Generally, the program reduces 
more quickly for smaller problem size, as would be expected. 

Storage . - Four terms are involved in considerations of storage for this program; 
they are related by the equation 

PROGRAM STORAGE + INITIAL EQUATION STORAGE 

+ Storage + g storage = C0NSTANT 
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0 100 200 300 ®0 500 600 

System size, number of equations 

Figure 15. - Comparison of reduction rate as function of problem 
size for complex and simple ladders. 


Since program storage is fixed when considering problem capability for a given computer, 
we have 


INITIAL EQUATION STORAGE + Xg TORAGJE + ^STORAGE = CONSTANT 
Now for large problems with relatively few inputs 


N e q = NUMBER OF EQUATIONS w NUMBER OF X's = N x 
and for a given class of problem we can write 

X STORAGE = K ’ N EQ 


Also, it is clear that 


INITIAL EQUATION STORAGE = K"N eQ 
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n - K*’ n N 

'■'STORAGE “ G 


Combining these results yields 

k ’ n eq + k ” n eq + k " tn g = CONSTANT 


Hence, 


N eQ + KN g = CONSTANT 


Here N„ Q is the original number of equations and Nq is the total number of G’s 
(original plus super G’s). This relation then should be useful in setting or changing 
array sizes for either (1) different computers or (2) problems where the relation of Nq 


*° n eq 


varies substantially. 


CONCLUSIONS 

Techniques have been established and a computer program written to allow the com- 
puter reduction of large arbitrary block diagrams. The techniques used involve symbolic 
manipulation of the defining equations by means of the FORMAC symbolic computer lan- 
guage. With the program, transfer functions can be determined for any output variable 
to any or all input variables. Since the primary objective of the work was to handle quite 
large and complex block diagrams, the solution is not formed as a single relationship for 
a desired transfer function. The solution, instead, is determined as a set of 
algebraic equations. Each of these equations, called super G’s, is defined in terms of 
the original system G's and the super G’s which precede it. 

The transfer functions are represented as a set of nested equations. This form has 
the advantage of being considerably more compact than the single relationship. Also, 
since the solution-set relations can be oulput from the computer as they are formed, the 
critical problem of computer storage is alleviated. Further, the nested form of the solu~ 
tion set allows straightforward numerical evaluation. 

The program is premised on the reduction of the block diagram by loops. When this 
is achieved, the solution is automatically scaled so that effective large powers of fre- 
quency do not occur. Also, the denominators are such that the possibility of zero- 
divides is small. 

Two forms of the solution set are possible at the option of the user. 

A program has been written to evaluate the solution set. This evaluation yields the 
real and imaginary parts and the phase angle and magnitude and log magnitude ratios of 
the transfer functions. 
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The advantages of the technique developed in the present study over previous tech- 
niques are as follows: 

1. As in the predecessor program, the arrangement of the input data in the form of 
arrays and the matrix manipulation for each frequency point are not required. (The 
symbolic solution set need be generated only once and then is reevaluated for each fre- 
quency. ) 

2. Unlike the predecessor program, the generation of a solution set instead of a 
single function allows the solution of considerably more complex problems, that is, sys- 
tems involving as many as 600 equations and a total of 1200 G’s. 

The user input to the reduction program consists of two basic pieces of information: 
(1) the algebraic equations representing the block diagram; and (2) the order string, 
which controls the manner in which the reduction is performed. The first part repre- 
sents a minimal effort on the part of the user. The determination of the order string 
does pose a small burden for the user. However, usually the consequence of a poor 
choice for the order string is a relatively small increase in the size of the solution set. 

The computing time required for the program varies with the size of the problem 
being solved. Typical rates are 0. 6 variable per second for small problems to an initial 
rate of 0.1 variable per second for the larger problems. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 5, 1973, 

501-24. 
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APPENDIX A 


SYMBOLS 

DENOM denominator of super G 

G transfer function for block in original block diagram 


g super 

function defining part of solution set 

i 


K 

constant 

N () 

number of ( ) 

S 

Laplace variable 

X 

signal or variable in a block diagram 

< X > 

circular frequency, rad/sec 

Subscripts: 


a,b, c, ♦ * . , 

index numbers 1, 2, . . . 



EQ 

equations 

i 

inputs 

OS 

order string 

V 

variables 
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APPENDIX B 


ORDER STRING 

The order string is a list of the block diagram variables. The order of the variables 
occurring in the list is the order in which they are to be reduced in the program. The 
variables are separated by commas. For example, X a> X^, X c indicates eliminate vari- 
able X , then X. , then X . Separation of a pair of variables by double commas indi- 
cates that a super G substitution is to be made at that point, 

The program logic is premised on reduction of the block diagram by loops. The 
loops of the block diagram are analogous to the cycles or circuits of graph theory (refs. 

4 and 5). That is, a loop is a progression of variables in a block diagram that closes on 
itself. Variables which split off (nodes) and leave the loop or those which enter the loop 
(through a summer) will be called communicating spurs, or simply spurs. 

The suggested best approach then will be to identify from the block diagram those 
loops having the least number of communicating spurs. The variables composing the 
loop, less the communicating spur variables, if possible, then form an order substring. 
The order of reduction of the variables forming the loop does not appear to be important. 
In the example of figure 16, the loop comprised of Xg, X 4 , X^q has only two spurs, 



Figure 16. - Block diagram used to show reduction of order string. 


namely Xg and Xg. Hence, an appropriate substring would be X 4 , X J0 or X 1Q , X 4< 
Elimination of this loop would leave a new loop with variables Xg, Xg, and X g with 
spurs Xg and X 2 giving a substring X g , Xg. After the reduction of each loop, a super 
G substitution should be made, that is, separation by double commas. Hence, the order 
string for this reduction might be 

X 10 , X 4 , , Xg, X 8 , , Xg, Xg, , X 2 , X 7 

where we have considered X g to be the output variable. 

While it is not always possible to get an end loop, that is, one with only two spurs. 
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the program appears to be quite forgiving; and even taking internal loops as shown m the 
text, gives rise to the same general forms. This will generate a larger number of super 
G's, however. 

Several checks can be made after an order string has been selected: 

(1) All variables in the problem must be represented in the order string except the 
output and the inputs (which, of course, are not eliminated in the reduction). 

(2) No variable can be represented more than once. 

(3) From algebraic theory, 


for the original system of equations. And, for the order string, 

(N ) = N - N. - 1 

v v'os V 1 


Hence 


<Vo S = n eq - 1 

This check is also made by the program. 

The choice of the order string for very large or very complex (interactive) diagrams 
may be difficult. Occasionally, for such cases the solution generated will be valid only 
over some limited frequency band. This is usually recognized by an instability occurring 
in the numerical solution above some frequency. In such cases, it may be necessary to 
change the order string to achieve a better solution form. 
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APPENDIX C 


PROGRAM IMPLEMENTATION 

The complete reduction and numerical evaluation of a block diagram is done in two 
separate programs. The reduction program is called REDUCE II and is written in FOR- 
MAC, and the evaluation program is called EVAL II and is in FORTRAN. FORMAC is 
an algebraic manipulation language. FORMAC stands for FOR mula MAn ipulation 
Compiler, It is an experimental language developed by IBM, and it runs on an IBM 
7094/7044 computer. Bond describes FORMAC in reference 3. There are other alge- 
braic languages which might be used in block diagram reduction. These languages are 
detailed by Sammet in reference 6. Two programs are necessary to maximize block- 
diagram-size capability. This approach separates the numerical evaluation from the 
reduction process. This separation is made since the numerical evaluations are per- 
formed more efficiently in FORTRAN than in FORMAC. The two-program approach 
also gives the user the ability to change component descriptions in EVAL n. 

The program listings and flow charts are given in appendix D. 


REDUCTION PROGRAM - REDUCE II 

REDUCE II is the main program in the block diagram reduction. Initially, the block 
diagram expressions, the inputs and the output, are read and counted. The reduction 
order string is examined for algebraic compatibility. That is, 

Number of variables in order string = Number of expressions - 1 

Since the user may solve for the ratio of an output to an input or an output to all the in- 
puts (see User’s Manual (appendix E) for details), the block diagram inputs not used in a 
specified transfer function are set to zero. 

The reduction proceeds by choosing a variable to eliminate, solving for that variable 
in one of the equations, and substituting the solution for that variable in the remaining 
equations. Subroutine INSTRN finds the variables in a reduction substring and puts their 
subscripts in common block STRING. Subroutine SEARCH picks the equations for elimi- 
nation and substitution in the reduction STRING. SEARCH obtains the variable to be 
eliminated from the information contained in common block STRING. The general rule 
is to take the simplest equation as the one to be eliminated. By simplest is meant the 
equation with the fewest variables. If the variable due for elimination is in more than 
two equations, some more selective rules are employed. These rules are embodied in 
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the following priority schedule, which has as its basis identification of the equations in 
the loop specified by the order substring (when possible): 

(1) If another variable in the reduction substring is contained in two or more equa- 
tions that contain the variable to be eliminated, pick the simplest equation of this new 
subset of equations to substitute into the other equations that contain the variable to be 
eliminated. 

(2) If there are no more variables in the reduction substring or if rule 1 is not satis- 
fied, look for other variables in common among the equations that contain the variable to 
be eliminated. If two or more equations contain a variable in common (other than the 
variable to be eliminated), pick die simplest of these as the equation to be eliminated. 

If neither rule 1 nor rule 2 can be followed, the simplest equation is picked as the one to 
be eliminated. Further, if only one (or no) equation contains the variable to be elimi- 
nated, an error message is written out and execution is stopped. With this method, the 
number of equations and variables is reduced by one at each step. The variables that 
appear in the transfer function are, of course, never eliminated. The final equation 
contains only the variables of the transfer function and the G-functions. 

In the reduction process, new combinations of G's (functions) are formed in sub- 
routine SUPERG. It should be noted that the program assumes that an expression is an 
equation. That is, if G(1)*X(1)+X(2) is the FORMAC expression, then G-jXj + X g = 0 is 
the implied equation. An expression prior to super G substitution might look like 

(1. +G(1)*G(2))*X(1)+G(3)*G(4)*X(5)-G(5)*X(6) 

The SUPERG subroutine would pick (l.+G(l)*G(2)) as a denominator and construct new 
coefficients such as G(6) and G(7), where 

G(6)=G(3)*G(4)/(1. +G(1)*G(2)) 


and 

G(7)=-G(5)/(l- +G(1)*G(2)) 

SUPERG would punch these new coefficients in FORTRAN in the following style: 

DENOM(l)=l. +G(1)*G(2) 

G(6)=G(3)*G(4) 

G ( 6 )=G (6 )/DENOM ( 1 ) 

G(7)=-G(5) 

G(7)=G(7)/DENOM(l) 
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Upon leaving subroutine SUPERG this expression would be X(1)+G(6)*X(5)+G(7)*X(6). 

The user has the option of selecting one of two forms for a super G. It is possible 
to always have a denominator of the form 

l-+VVVVV G f*VVV G m' . . 

This can be done by picking a coefficient of some variable and adding 1.0 to it. This new 
term can be divided into all the coefficients of the variables in an expression. This pro- 
cess tends to generate more super G's than when the super G’s are allowed to form 
naturally. By punching the word ARTIFICIAL on the first data card to REDUCE II this 
form of super G will be produced whenever SUP ERG finds a coefficient with two or more 
terms. If the user does not specify ARTIFICIAL or NATURAL (1.0 is not added to coef- 
ficients), an error message is printed out and execution stops. 

The output of the REDUCE II program is punched on cards and printed on paper. 

The cards plus the user’s data supply all the information for the EVAL II program. See 
the user’s manual for further information of card input and output. Since the program 
should be as easy to use and as reliable as possible, the program input is given in alge- 
braic terms for the most part. The user does not have to give a count of how many vari- 
ables, expressions, inputs, and functions are to follow, as was required in the predeces- 
sor of this program. The programs interpret the user’s data for this information. The 
only nonalgebraic piece of data is the first card, which specifies what form the user de- 
sires the super G’s to have. An attempt has been made to make the programs efficient 
in terms of execution time. The subroutine COUNTX was written to be faster than the 
corresponding FORMAC routines that could have been used. Auxiliary array sizes have 
been kept as small as possible to maximize block-diagram-size capability. It is possible 
to have a block diagram that is so large that it would be stored on an auxiliary storage 
device (tape or disk file). If this does happen, there will be a substantial increase in 
execution time because moving information to and from an auxiliary device is inherently 
slow. REDUCE II does not call on SIN, COS, TANH, ATAN, ALOG, DUMP, PDUMP, 
FMCDIF, EXPDMP, or FMCDMP. Using the dummy subroutines SIN, COS, TANH, 
ATAN, ALOG, DUMP, PDUMP, FMCDIF, EXPDMP, and FMCDMP allows us to use the 
storage space normally assigned to these routines. 

The user may, for any particular job, change any of the dimensions of the program. 
In particular, if more inputs are desired, the array INPUT should be lengthened. If 
more equations or variables are needed, the arrays EQ and X should be changed. Also 
the value of MAX(l) must be changed. If some other number of functions (G’s) is to be 
used, the G array and MAX(2) should be changed accordingly. 
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EVALUATION PROGRAM - EVAL II 


The EVAL II program evaluates the cards REDUCE II punched out. The user in- 
serts some of these cards directly into EVAL II as FORTRAN coding. The remainder is 
numerical data to identify transfer functions. The program covers a frequency band of 
FSTART to FEND in steps of DELTAF. The output of EVAL II consists of a tabulation 
of the transfer functions for each frequency point. Also a set of CALCOMP plots for the 
frequency band can be obtained. The written information gives a transfer function num- 
ber, the frequency point, the real and imaginary parts of the transfer function, the abso- 
lute value and phase angle of the transfer function, and the magnitude of the transfer 
function in decibels, based on dB = 20 log^Q (absolute value). The plots are of magni- 
tude against frequency and phase angle against frequency. These plots can be made on a 
linear or semi logarithmic set of axes. The user’s manual (appendix E) gives more in- 
formation on the actual use of EVAL II. 


SUBPROGRAMS 


REDUCE II uses the following subprograms that are not supplied with either FOR- 
TRAN, FORMAC, or this report. These subprograms must be supplied by the user: 


Function 
ALS(N,X) 
IARS(N, X) 
LGR(N, X) 
AND(Xl, X2) 
OR(Xl, X2) 


Description 

Accumulator left shift of X by N binary places 
Accumulator right shift of X by N binary places 
Logical right shift of X by N binary places 
Logical intersection of XI and X2 
Logical union of XI and X2 


EVAL II uses the following subprograms that are part of the CALCOMP plotting 
package. These programs must be user supplied since CALCOMP has proprietary 
rights on their use. They are LINE, SCALE, AXIS, SYMBOL, NUMBER, and PLOT. 
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APPENDIX D 


DESCRIPTION, FLOW CHARTS, AND LISTINGS FOR COMPUTER PROGRAMS 

DESCRIPTION 

Reduction Program - REDUCE II 

The main program REDUCE II reads the input data, interprets data, and does the 
actual block diagram reduction. It prints and punches out the transfer function gener- 
ated. 

Subroutine COUNT counts the number of variables in an expression. It returns the 
number of variables and/or the variables to the calling program. 

Subroutine SEARCH finds expressions containing the variable to be eliminated. It 
constructs a list of the variable to be eliminated, the expression to be eliminated, and 
the equations to be substituted into. 

Subroutine SUPERG creates the super G’s or combinations of functions formed in 
the reduction process. 

Subroutine SOUT outputs FORMAC expressions. This output consists of a printed 
statement, a punched card, or both. 

Subroutine INSTRN reads reduction substrings and provides a count of the number 
of variables and the variables in the substring. 

DELETE, REMOVE, FC, and NBR are used to remove $ from FORMAC expres- 
sions and to convert floating-point exponents to fixed-point exponents. 

Subroutine COS with entry points SIN, TANH, ATAN, ALOG, DUMP, PDUMP, 
FMCDIF, EXPDMP, and FMCDMP is a dummy subroutine that is not called. It pro- 
vides entry-point names that are the same as subprograms referenced by FORMAC but 
are not executed by the programs of this report. This prevents the equivalent subpro- 
grams from being loaded and conserves storage space. 


Evaluation Program - EVAL II 

The main program EVAL II reads and writes all the numerical input and output. It 
calls the subprograms that evaluate the super G's and the transfer functions. 

The function subprogram EXPFX, SINFX, COSFX, SINHFX, COSHFX, SQRTFX 
evaluates the supplied system functions. 

The G's, super G’s, and transfer function cards (NO=, DO=) are inserted in sub- 
routine COEFF. The actual evaluation of the transfer function occurs in COEFF. 
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SYMARG A 

ATOMIC XI631 I f GI1202) A 

COMMON /TRANS/ ART IF 

DIMENSION INI70I, OUTI17I, NEQUOOJ, EQI601I, INPUTI5),. NXARRYI30) A 
If MINXGI 2) t GORX( 2) t MAX ( 2 1 1 LSAVE I 50 I A 

DIMENSION IARTIFI 3) .NATURI3I 

LOGICAL LOCATE, TEST, FINISHt DONEt ART IF A 

INTEGER OUT 

REAL ISfITtNUM A 

EQUIVALENCE ( INI 1 ) , OUT 1 1 ) , N£Q( III, 1 MINXGI 1 ) » Ml t I HI NXS I 21 ,MG I A 

EQUIVALENCE I PLUS1, ARTIF I 

DATA IARTIFI 1), NAT UR ( 1 ) /15HART IF I Cl AL FORM, 13HNATURAL FORM / 

DATA MAX/630,1200/ A 

DATA DR/6HAOOOOO/ A 

DATA XX/6HKOOOOO/ A 

DATA MINXG,G0RX/2*0,1HX,1HG/ A 

C**« '******************************************************************** A 
C M IS THE NUMBER OF VARIABLES A 

C Ml IS THE NUMBER OF EQUATIONS A 

C IK IS THE NUMBER OF INPUTS A 

C MG IS THE NUMBER OF FUNCTIONS AT INPUT A 

C **************** 4 *** *************** ******** ******* ************* ******** A 
1 CONTINUE A 


PLUS1=3. 

REA0(5,42> (INIII, 1=1,31 
00 100 1=1,3 

IF! INI II. EQ. IARTIFI III PLUS1*PLUS1*1. 

IFC INI I I.EQ.NATURt I ) I PLUS1=PLUS1-1. 

103 CONTINUE 

IFIPLUSl.NE.O. . AND.PLUS1.NE.6. I GO TO 300 
C******************** *************************************************** A 


C READ THE BLOCK DIAGRAM PARTS A 

£*« ********************** ************************* ********************** A 
NX = 0 A 

Ml = l A 

IK=0 A 

MG=0 A 

M=0 A 

DO 19 J= 1 , A A 

NN=0 A 

PREV=BL A 

2 K=1 A 

3 K1 = 14*1 K-l I +1 A 

K2=K1+13 A 

READ (5,421 I INIL I ,L=K1,K2) A 

GO TO (4,9, 13,171, J A 

C **************** ******************************************************* A 

C INTERPRET THE BLOCK DIAGRAM INPUTS A 

C ****************************************** 4 *******♦*♦♦**♦♦♦**♦* **♦***♦♦ A 

4 CONTINUE A 

DO 5 1*1,80 A 

F*FC ( IN( KI I , IJ A 

IF <F*EQ.DM IK-1K+1 A 

IF (F#EO*DR#AND,PREV«EO*DR) GO TO 6 A 

PREV*F A 

5 CONTINUE A 

K*KU A 

GO TO 3 A 


1 

2 

3 

4 

5 

6 
7 


8 

9 

18 

11 

12 

13 

14 

15 
IS 
17 
21 


10 

19 

20 
22 

23 

24 
24A 
24B 

25 

26 
27 
20 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 
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JJ = D 
IKMK-i 

IF < IK.GT.5) GO TO 35 
DO 7 JJJ*1, IK 

LET I NPUT ( JJJ > = ALGCON IN<1),JJ 

7 CONTINUE 
WRITE (6,54) 

DO 8 L = 1 , IK 
Q=0.0 

LET 0 = 8CDC0N INPUT ( L ) , OUT ,17 

8 WRITE (6,51) OUT( 2 I 

GO TO 19 

C *********************************************************************** 

C INTERPRET THE DESIRED BLOCK DIAGRAM OUTPUT 

C *********************************************************************** 

9 CONTINUE 
JQ*0 

DO 10 1*1,80 

10 IF (FC( IN(K1 ),1 J.EQ.DR ) JQ=JQ*1 

JJ*0 

LET NUM = ALGCON IN(i),JJ 
IF ( JQ. EQ* 2 ) GO TO XI 
LET DEN * ALGCON IN(1),JJ 

11 CONTINUE 
Q=0.0 

LET Q * BCDCON NUM,OUT,17 
IF (JQ.EQ.3) GO TO 12 
WRITE (6,46) OUT! 2 ) 

WRITE (6,49) 

GO TO 19 

12 CONTINUE 

WRITE (6,47) OUT( 2 ) 

Q=0.D 

LET Q * BCDCON DEN, OUT, 17 
WRITE (6,48) OUT ( 2 I 
WRITE (6,49) 

GO TO 19 

C * **** ** ** *** * * ***** * ** ** **** ********* ** * **** * **** ** * ***** * * ** ******* * * * 
C INTERPRET THE BLOCK DIAGRAM EXPRESSIONS 

C ******************************************************* **************** 

13 CONTINUE 
Q*l* 

DONE** FALSE ♦ 

DO 14 1=1,80 
F*FC ( IN( Kl> , I) 

IF (F*EO.OR) Q=0. 

IF (F.EG.DR.ANO.PREV.EQ.DR) DONE*. TRUE. 

PREV=F 

14 CONTINUE 

WRITE (6,53) Ml , { IN(KK),KK=K1,K2) 

K=K+1 

IF (K.GT.5) GO TO 38 
IF (Q.NE.O* ) GO TO 3 
DO 16 MM=1,2 

CALL COUNT ( ! N( 1 ) , GORX(MM) ,NXS,NXARRY,0) 

IF (NXS.EQ.O) GO TO 16 
DO 15 NN=1 , NXS 

15 MI NXG( MM) =MAX0( MI NXG( MM ) ,NXARRY( NN) I 
IF (MINXG(MH).LE.MAX(MM) ) GO TO 16 
GO TO (33,32) , MM 


A 45 
A 46 
A 47 
A 48 
A 49 
A 53 
A 51 
A 52 
A 53 
A 54 
A 55 
A 56 
A 57 
A 58 
A 59 
A 60 
A 61 
A 62 
A 63 
A 64 
A 65 
A 66 
A 67 
A 63 
A 69 
A 70 
A 71 
A 72 
A 73 
A 74 
A 75 
A 76 
A 77 
A 78 
A 79 
A 80 
A 81 
A 82 
A 83 
A 64 
A 85 
A 86 
A 87 
A 88 
A 89 
A 90 
A 90 
A 91 
A 92 
A 93 
A 94 
A 95 
A 96 
A 97 
A 98 
A 99 
A 100 
A 101 
A 102 
A 103 
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16 CONTINUE 

IF (Ml.GT.MAXIin GO TO 34 
JJ = D 

LET EQI Ml )* ALGCON 1N(1)»JJ 
LET EQ( Ml )®EXPAND EG(Ml) 

M1 = M1U 

IF (.NOT. DONE) GO TO 2 
Ml=Ml-l 

C***** ************************************************************ ****** 
C INTERPRET AND WRITE THE BLOCK REDUCTION ORDER 

L***** ************ ******* ******** *************************************** 


17 NN=NNU 

IF (NN.EQ.l) WRITE (6,44) 

WRITE (6,53) NN»( IN(L),L=K1,K2) 

DO 18 1*1, BO 
CHAR=FC< IN(K1), I) 

IF (CHAR. EG. XX) NX=NX+1 

10 IF (CHAR. EO. OR) GO TO 19 

GO TO 2 

19 CONTINUE 

IF (NX.NE.M1-1) GO TO 37 

C *********************************************************************** 

C BACKSPACE THE INPUT UNIT 

C ********************************************* ************************** 


DO 20 L=1»NN 
20 BACKSPACE 5 

WRITE (6,45) M, Ml , MG 
WRITE 16,39) 

Q ***** *************** *************************************************** 

INITIALIZE AND set UP EQUATIONS FOR SOLUTION 


IF ( JQ.EQ.2 I GO TO 23 

REPLACE THE tNPUTS BY ZERO IF NOT DEN OR NUM 


DO 22 L=1»M1 
DO 21 J=1,IK 

LET TEST = MATCH ID, NUM, INPUT( J ) 

IF (TEST) GO TO 21 

LET TEST = MATCH ID, DEN, INPUT! J) 

IF (TEST) GO TO 21 

LET EQ(L) = SU8ST EQ(L ) , ( INPUT ( J ) , 0 1 

21 CONTINUE 

22 CONTINUE 

c *********************************************************************** 

■>•» rnwTiNHF 


C REDUCE THE EQUATIONS 

c*********************************************************************** 

ot rnwriMiiF 



C GET A VARIABLE TO ELIMINATE 

C 

25 CONTINUE 


A 104 
A 105 
A 106 
A 10? 

A 109 
A 109 
A 110 
A 111 
A 112 
A 113 
A 114 
A 115 
A 115 
A 117 
A 118 
A 119 
A 120 
A 121 
A 122 
A 123 
A 124 
A 125 
A 126 
A 127 
A 120 
A 129 
A 130 
A 131 
A 132 
A 133 
A 134 
A 135 
A 136 
A 137 
A 130 
A 139 
A 140 
A 141 
A 142 
A 143 
A 144 
A 145 
A 145 
A 147 
A 143 
A 149 
A 153 
A 151 
A 152 
A 153 
A 154 
A 155 
A 156 
A 157 
A 158 
A 159 
A 163 
A 161 
A 162 
A 163 
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CALL SEARCH (X, EQ, M, Ml ,NEQ, DONEI 
11 = 1 

J=NEQU1) 

IF (J.EQoOl GO TO 27 
C 

C GET AN EQUATION CONTAINING THE VARIABLE 

C 

11 = 11+1 
L = NEQ( 1 1 ) 

C 

C SOLVE THE EQUATION 

C 

LET RR * EQ< L I 
LET EQ( L)=0« 

LET IT * COEFF RR,XU> 

LET RR * EXPANO RR-IT*X(J) 

C 

C ELIMINATE THE VARIABLE IN THE REMAINING EQUATIONS 

C 

26 11=11+1 

L=NEQ( II) 

If (L.EO.O) GO TO 27 

LK=LK+1 

LS AVEI LK >=L 

LET IS = COEFF EQ(L>»X(J) 

LET EQID = EXPAND EQIL >-IS*X(J> 

LET EO(L) = EXPAND IS*RR-IT*EQ< L I 
GO TO 26 

27 CONTINUE 
ERASE IS, IT 

IF (.NOT.DDNEI GO TO 25 
IFIFINISHI GO TO 280 
DO 28 MM*l,LK 
NN=LK-MM+1 
L=LSAVE«NN» 

C ***** ************ ****************************************************** 

C CREATE SOME SUPER G'S 

C*^****************^*********^^t*************^m*m*******ti***^**^**« 

CALL SUPERG <EQ (1 ) , NG.G, X , MAXI 2>, M I 

28 CONTINUE 
GO TO 2* 

C 

C 

280 L=NEQ( 1 1-1 1 

C SOLVE FOR NUMERATOR AND DENOMINATOR 

C **************** ******************************************************* 

LET R * EQ( L I 
LET RD = COEFF R, NUM 
LET EQ(L)=EQ(L)-RD*NUM 
LET EQIL > =EXPAND EQILI 
LET RD = EXPANO -1*RD 
PUNCH 61 

CALL SOUT ( RD, 2, 1H ) 

CALL COUNTX < NUM, IHX, NXS, NX ARR V,0) 

NNUM=NXARRV( 1 ) 

IF (JQ.EQ,3) IK=1 
DO 30 1=1, IK 
IF (JQ.EQ.3) GO TO 29 
LET OEN=INPUT( I » 


A 16* 
A 165 
A 166 
A 1ST 
A 169 
A 169 
A 175 
A 171 
A 172 
A 173 
A 17* 


A L77 ' 
A 179 
A 179 
A 183 
A 191 
A 192 
A 183 
A IB* 

A 185 
A 186 
A 187 
A 188 
A 199 
A 195 
A 191 
A 192 
A 193 
A 19* 

A 195 
A 196 
A 197 
A 198 
A 199 
A 205 
A 201 
A 202 

A 20* 

A 205 

A 207 
A 203 
A 209 
A 210 
A 211 
A 212 

A 21* 

A 215 

A 217 
A 218 
A 225 
A 222 
A 223 
A 22* 
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29 CONTINUE A 225 

CALL COUNTX { DEN, 1HX, NXS, NXARRY, 0) A 226 

OUTm=NXARRY(l> A 227 

LET RN = COEFF R.DEN A 229 

LET EQ(L>=EQ(L)-RN*DEN A 233 

LET EQ(L >=EXPAND EQ(L) 

C ******* ********* ****** ************************************************* A 232 
/C WRITE NUMERATOR AND DENOMINATOR A 233 

' C*********************************************************************** A 234 
WRITE (6, A3) NNUM, OUT ( 1 1 A 235 

WRITE (6,50) A 236 

PUNCH 60, I A 237 

CALL S0UT(RN,3,1H ) 

WRITE (6,52) A 239 

CALL SOUT (RD,1,IH ) A 2 A3 

30 CONTINUE A 241 

CALL COMOUT A 242 

PUNCH Al, NNUM, IK, (OUT! I)»I=l,IK) 

ERASE RN,RD A 2 A3 

C *********************************************************************** A 244 

C CHECK ON THE COMPLETEDNESS OF THE SOLUTION A 245 

C ***** ****************************************************************** A 246 
LET R=0. A 247 

DO 31 L=1»H1 A 248 

LET TEST= MATCH 10,EO(L),R A 249 

IF (TEST) GO TO 31 A 250 

WRITE (6,62) L,L A 251 

CALL SOUT ( EQ( L ), 1, 1H > 

31 CONTINUE A 253 

GO TO 1 A 254 

£******************«******************************************** ******** a 255 

C WRITE ERROR MESSAGES A 256 

C **************** ******************************************************* A 257 

32 CONTINUE A 258 

WRITE (6,57) MAX( 2) A 259 

GO TO 36 A 260 

33 WRITE (6,55) MAX(l) A 261 

GO TO 36 A 262 

34 WRITE (6,56) MAX(l) A 263 

GO TO 36 A 264 

35 WRITE (6,58) A 265 

36 STOP A 266 

37 WRITE (6,40) A 267 

GO TO 36 A 268 

38 CONTINUE A 269 

WRITE (6,59) A 270 

GO TO 36 A 271 

380 WRITE(6,63) 

STOP 

C A 272 

C A 273 

39 FORMAT (//IOX.33H* * * OUTPUT FROM REDUCE II * * M A 274 

40 FORMAT ( 10X , 114H**ERR0R** NUMBER OF VARIABLES IN REDUCTION STRING A 275 

IIS NOT CONSISTENT WITH NUMBER OF EXPRESSIONS — EXECUTION STOPPED ) A 276 

41 FORMAT (15 ) A 277 

42 FORMAT (13A6,A2> A 273 

43 FORMAT ( 1HK , 10X ,20HTRANSFER FUNCTION X( , T3,4H) /X ( , I 3, 2H ) ) A 279 

44 FORMAT (23HKTHE REDUCTION ORDER IS,//) A 283 

45 FORMAT (34H0THE LARGEST VARIABLE SUBSCRIPT IS,I5,/,29H THE NUMBER A 2BI 

10F EXPRESSIONS IS,I5,/,38H THE LARGEST FUNCTION (G) SUBSCRIPT IS, I A 282 
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/ 



24) 


A 

283 

*6 

FORMAT 

I21HKTHE SOLUTION IS FDR ,A6,14H TO THE INPUTS! 

A 

284 

47 

FORMAT 

(21HKTHE SOLUTION IS FOR ,A6,2HT0> 

A 

285 

48 

FORMAT 

< LH*,29X,A6) 

A 

266 

49 

FORMAT 

( 30HKTHE BLOCK DIAGRAM EXPRESSIONS,//) 

A 

287 

53 

FORMAT 

C1HX,15X,13HTHE NUMERATOR) 

A 

288 

51 

FORMAT 

( 5X, 9A6) 

A 

289 

52 

FORMAT 

(IHK,15X,15HTHE DENOMINATOR ) 

A 

293 

53 

FORMAT 

<3X, 14, 2X, 13A6, A2 ) 

A 

291 

54 

FORMAT 

(25HITHE BLOCK DIAGRAM INPUTS,//) 

A 

292 

55 

FORMAT 

(20H **ERROR** MORE THAN,I4,10H VARIABLES) 

A 

293 

56 

FORMAT 

(20H **ERROR** MORE THAN,I4,10H EQUATIONS) 

A 

294 

57 

FORMAT 

<20H **ERROR** MORE THAN,I4,10H FUNCTIONS) 

A 

295 

58 

FORMAT 

<43H **ERROR** MORE THAN 5 BLOCK DIAGRAM INPUTS) 

A 

296 

59 

FORMAT 

(59H **ERROR** MORE THAN 5 CAROS FOR A 8LOCK DIAGRAM EXPRES 

A 

297 


IS I ON > 


A 

298 

60 

FORMAT 

(6X,3HN0(,I2,2H)*> 

A 

299 

61 

FORMAT 

(6X,3HD0=) 

A 

303 

62 

FORMAT 

(//, 10X.69H***** THE SOLUTION FOR THIS PROBLEM MAY BE IN ER 

A 

301 


1ROR BECAUSE EQUATION , I4.31H WAS NOT ELIMINATED, EQUATION , I4,7HI 

A 

302 


2S*%***> 


A 

303 

63 

FORMAT 

C10X, 40HINVALID OPTION FOR THE SUPER G FORM CARD) 




END 


A 

304- 
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subroutine countx < eq»gorx,nxs»nxarry, itype i B 1 

ENTRY COUNT G ( EQ,NXS»GORX, ITYPE ) 8 2 

SYMARG EG 83 

DIMENSION NXARRY41I, MTEMP4 4 >» OUT415), TEMP44I, OTEMP414) B 4 

DIMENSION FC241I 85 

EQUIVALENCE (MTEMP, TEMPI 8 5 

EQUIVALENCE (OUT( 2 1 tOTEMPI l II B 7 

LOGICAL NBR, DONE B B 

LOGICAL ZERO 6 9 

DATA AST, DL,RP/6H*00000»6H*00000,6H) 00000/ B 13 

NXS=0 B 11 

ICODE^O B 12 

DONE®. TRUE* 8 13 

IU®84 B 14 

QJ=0, B 15 

1 CONTINUE . B IS 

LET QJ=BCOCON EQ»0UT»15 B 17 

GO TO 4 '' 8 18 

ENTRY COUNT 4 ARRAY, GORX,NXS,NXARRY, I TYPE I B 19 

DIMENSION ARRAY <11 B 23 

IC0D£=1 B 21 

DONE*, TRUE, B 22 

NXS=Q B 23 

IU»80 B 24 

K=0 8 25 

2 K=K>1 B 26 

K1 = 4 K.-1 ) *14+1 B 27 

K2=K1+13 B 28 

L=0 B 29 

DD 3 I*Kl»K2 - B 33 

L=L*1 B 31 

3 OTEMP4L)=ARRAY4 II B 32 

4 CONTINUE B 33 

IF 4 .NOT *DONE ) GO TO 5 B 34 

AX=FC14G0RX,0) B 35 

PREV1-AX B 36 

5 CONTINUE B 37 

J=0 B 38 

6 J=J*1 B 39 

IF I, NOT. DONE) GO TO 8 B 40 

EC2=FC4 OTEMP, J) B 41 

IF 4 ITYPE. NE.l) GO TO 7 B 42 

IF 4 PREV1.EQ. AST.AND.FC2.EQ* AST ) NXS=NXS+1 B 43 

PREV1=FC2 B 44 

7 CONTINUE B 45 

IF 4FC2.NE.AXI GO TO 14 B 45 

DONE*. FALSE. . B 47 

NXS=NXS+1 " B 48 

KTR=0 B 49 

J=J*2 B 53 

IF (J.GT.IUI GO TO 15 B 51 

8 CONTINUE B 52 

DO 10 J2=J» XU B 53 

J3=J2 B 54 

FC2=FC40TEMP,J2> B 55 

IF 4FC2.EQ.RP) GO TO 9 B 56 

IF 4.NOT.N8R4FC2, III GO TO 10 B 87 

KTR=KTRU 0 59 
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TEMP(KTRI=ARS<30,FC2> 

IF (J2.EQ.IU) GO TO 15 
GO TO 10 

9 CONTINUE 
DONE=.TRUE. 

GO TO 11 

10 CONTINUE 

11 CONTINUE 
J=J3 

IF I. NOT. DONE) GO TO 15 
IF (ITYPE.EQ.il GO TO 14 
NXARRY ( NXS) =0 
ITEN=1 

DO 12 1 1 = 1, KTR 
KTR1=KTR-II+1 

NXARRY(NXS>=ITEN*MTEMP(KTR1)+NXARRY(NXS) 

12 ITEN=ITEN*10 

IF ( NXS* EQ. 1 1 GO TO 14 

NXSl=NXS-l 

DO 13 11=1, NXS1 

IF I NXARRYI NXS) .NE.NXARRY 1111) GO TO 13 

NXS=NXS-1 

GO TO 14 

13 CONTINUE 

14 CONTINUE 
DONE*. TRUE. 

IF (FC2.EQ.DL) GO TO 16 
IF (J.LT.IJ) GO TO 6 

15 CONTINUE 

IF (1C0DE.EQ.1I GO TO 2 
IF (QJ.NE.3.1 GO TO 1 

16 CONTINUE 
RETURN 
END 


B 59 
B 63 
B 61 
B 62 
B 63 
B 64 
B 65 
B 66 
B 67 
B 6B 
B 69 
B 73 
B 71 
B 72 
B 73 
B 74 
B 75 
B 75 
8 77 

B 78 
B 79 
B 83 
B 81 
8 82 
8 83 

8 84 

B 85 
8 86 
B 87 
8 88 
8 89 

8 93 

8 91 

B 92- 
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Subroutine SEARCH 
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SUBRUUdNE SEARCH ( X,EQ,MX,MEQ,NEQ,DQNE) 
DIMENSION NEO(l), NARRAY ( 30 ) * ESM 1) 

COMMON /STRING/ IA(200),M 
SYMARG X»EQ 
ATOMIC XU) 

LOGICAL DONE » Q 
DATA 11/0/ 

DONE*. FALSE. 

II = l I ♦ l 

NEQ( 1 )*IA( l 1 1 

L=0 

M I N-999 
DO 2 1*1, MEO 

CALL COUNTX ( EQ ( 1) , 1HX, NXS, N ARRAY, 0) 

IF INXS.EQ.O) GO TO 2 
DO 1 J=l» NXS 
K=NAKRAY( J) 

IF IK.NE.NEQU) I GO TO 1 
L=L + 1 

NEQ( L+2 ) = I 

IF (NXS. GE. MIN) GO TO 2 

MI N=NXS 

IMIN*I 

GO TO 2 fi'. 

CONTINUE 

CONTINUE 

IMSAVE= IMIN 

IF (L-2) 3, A, 9 

ERROR ONLY ONE EQUATION WITH VAR X(NEQU)) 
WRITE (6,17) NEQ(l) 

STOP 

NEQ( 2 ) = IMIN 
L1=L+1 
LH=L+2 
DO 6 1=3, LI 

IF ( NEQ( I ) . NE« IMIN) GO TO 6 
11 = 1+1 

DO 5 J=U,L2 
NEQ( J-1)=NEQ( J) 

GO TO 7 
CONTINUE 
CONTINUE 
NEQ(L2)=0 

IF ( Il.NE.M) GO TO 8 
11=0 

DONE*. TRUE. 

CONTINUE 

RETURN 

IF (II.EQ.M) GO TO 13 

11=11+1 

DO 12 1=11, M 

IS = 0 

MI N=999 

00 11 J=l,L 

K=NEQ( J+2) 

CALL COUNTX < EQ ( K ) , 1HX.NXS, NARR AY, 0 ) 

DO 10 N=i,NXS 

IF (N ARRAY! N).NE.IA(I)) GO TO 10 

0 


c 

1 

c 

2 

c 

3 

c 

4 

c 

5 

c 

6 

c 

T 

c 

8 

c 

9 

c 

10 

c 

11 

c 

12 

c 

13 

c 

n 

G 

15 

f . . ■ c 

IS 

c 

17 

c 

IB 

G 

19 

c 

20 

, ■ , • G 

21 1 

G 

22 

C 

23 

C 

24 

, ' ' ,c 

25 

c 

25 

c 

27 

.c 

28 

c 

29 

* , c 

30 

c 

31 

c 

32 

c 

33 

. c 

34 

, c 

35 

c 

35 

* ’ c 

37 

c 

38 

c 

39 

c 

40 

c 

41 

c 

42 

c 

43 

c 

44 

c 

45 

c 

45 

c 

47 

c 

48 

G 

49 

c 

50 

c 

5i 

c 

52 

c 

53 

c 

54 

c 

55 

c 

55 

c 

57 

c 

58 
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IS=1$+1 

c 

59 


IF INXS.GE.MIN) GO TO 11 

c 

60 


Ml N=NXS 

c 

61 


IMI N=K 

c 

62 


GO TO 11 

c 

63 

10 

CONTINUE 

c 

6* 

li 

CONTINUE 

c 

65 


IF < IS.LT.2 1 GO TO 12 

c 

66 


GO TO 4 

c 

67 

12 

CONTINUE 

c 

68 

13 

L1»L-1 

c 

69 


00 16 1=1, LI 

c 

73 


K=NEQ( I +21 

c 

71 


CALL COUNTX ( EQIK) , 1HX»NXS »NARRAY, 0 I 

c 

72 


J=I+1 

c 

73 


DO 15 N* J»L 

c 

7* 


Kl=NEQIN+2> 

c 

75 


DO 14 Nl*l, NXS 

c 

76 


N2=NARRAY(N1 ) 

c 

77 


IF ( N2» EQ.NEQI 1 I ) GO TO 14 

c 

79 


LET Q= FIND EQI K1 ) , APP, ALL, < X< N2H 

c 

79 


IF UNOT.Q) GO TO 14 

c 

83 


CALL COUNTG ( EQIKl ) ,NX1, IHX, 1 ) 

c 

81 


IMI N*K 

c 

82 


IF INX1.LT.NXSI 1 MIN=K1 

c 

B3 


GO TO 4 

c 

84 

1* 

CONTINUE 

c 

85 

15 

continue 

c 

86 

16 

continue 

c 

87 


IMIN=IMSAVE 

c 

68 


GO TO 4 

c 

89 

C 


c 

90 

17 

FORMAT (10X.36H**** ONLY ONE EQUATION CONTAINING X(,I4,6H) ♦***» 

c 

91 


END 

c 

92 
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Subroutine SUP ERG 



56 
















SUBROUTINE SUPERS ( EQ,MG,G,X,MAXG,MX) 
SYMARG EQ,X,G 
ATOMIC IllliGIll 
DIMENSION N ARRAY! 50 ) 

DIMENSION I SAMEI50) 

COMMON /TRANS/ ART IF 
LOGICAL TROTH 
LOGICAL FIRST 
LOGICAL PRESG 
LOGICAL YES 
LOGICAL ART IF 
DATA FIRST/. TRUE./ 

DATA MD/0/ 

IF (MG.GT.MAXG) RETURN 

LET EQSAVE*EQ 

MAX*~1 

COUNT THE NUMBER OF X'S IN EQUATION EQ 
CALL COUNTX ( EQ»IHX,NXS« NARRAY, 0) 

IF (NXS.EQ.O) GO TO 19 
DO 5 K=l,NXS 
l =NARRAY( K) 

LET RR=CO£FF EQ,Xf I ) 

LET IQ=CENSUS RR,TERM 

LET YES* FIND RR,APP,ALL,( 1.1 

IF (YES) GO TO 1 

LET YES* FIND RR, APP, ALL, 4 -1. I 

IF (YES) GO TO 1 

I Q*“IQ 

CONTINUE 

IF (IABSUQ).EQ.l) GO TO 5 
IF (MAX.GT.O.AND.IQ.GT.O) GO TO 2 

IF (MAX.LT.O.AND. IQ.GT.O) GO TO 3 

IF (MAX.LT.O.AND. IQ. LT.O) GO TO 4 

GO TO 5 

IF ( IQ.LT.MAX) GO TO 5 

HAX*IQ 

L-I 

GO TO 5 

IF I IQ.LT.MAX I GO TO 3 
CONTINUE 

IF (MAX.NE.-l) GO TO 7 
DO 6 K=1,NXS 
IsNARRAYIK) 

LET RR-COEFF EQ,X(I) 

CALL COUNTG ( RR.NGS, 1HG, 1 ) 

IF 4 NGS.LT . 3 ) GO TO 6 
FORM A NEM SUPER G 
MG*MG*1 

IF (MG.GT.MAXG) GO TO 20 
L=1 

LET EQ * EQ ~ RR * X(L) 

LET £Q*EXPAN0 EQ+ G(MG)*X4L) 

WRITE (6,24) 

WRITE (6,21) MG 
PUNCH 21, MG 
CALL SOUT( RR, 3, 1H+) 

CONTINUE 
GO TO 19 


D 1 
D 2 
D 3 
0 4 

D 5 

D 6 
D 7 
D 8 
D 9 

0 10 
0 11 
D 12 
D 13 
D 14 
D 15 
D 16 
0 17 

D 18 
D 19 
D 20 
D 21 
D 22 
0 23 

D 24 
D 25 
D 26 
D 27 
D 28 
D 29 
0 30 

D 31 
D 32 
D 33 
D 34 
0 35 

0 36 

D 37 
D 38 
D 39 
D 40 
D 41 
D 42 
D 43 
D 44 
D 45 
D 46 
D 47 
D 48 
D 49 
D 50 
D 51 
0 52 


D 54 
D 55 
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C DIVISION SEARCH AND CONSTRUCTION Of SUPER G»S 

7 CONTINUE 
PRESG*. FALSE. 

LET RR*COEFF EQ,X(L) 

LET EQ=EQ-RR*X(U 
LET DD= 1 • 

IF (.NOT.ARTIF.OR.MAX.GT.OI GO TO 8 
LET SS=RR 
LET RR=RR*1 * 

PRESG=.TRUE. 

8 CONTINUE 
M0=MD*-1 
WRITE (6,24) 

PUNCH 22, MD 
WRITE (6,22) MO 
CALL S0UT(RR,3, 1H>> 

IF ( .NOT .PRESG) GO TO 9 
MG=MG+1 

IF (MG.GT.MAXG) GO TO 20 
WR1TE(6,24) 

WRITE(6,21) MG 
CALL SOUT t SS , 1 , 1H+) 

LET DD=G(MG> 

PUNCH 27 , MG , MD»MD 
WRITE ( 6, 23) MG, MG, MD 

9 CONTINUE 
KTR=0 

PRESG*. FALSE. 

LET EQ-EXPAND EG+DD*X(L) 

KTR1=KTR 
DO 18 1*1 *NXS 
N=NARRAV( 1) 

IF ( N.EQ.L) GO TO 18 
IF (KTR.EQ.O) GO TO 11 
DO 10 J=1,XTR 

10 IF (N.EO.IABSt I SAME ( J ) ) ) GO TO 18 

11 KTR1-KTR1+1 
ISAME< KTR1) =N 

LET SS*COEFF EQ,X(N) 

DD 15 J*1,NXS 
K=NARRAY I J) 

IF (R.EQ.L.OR.K.EQ.N) GO TO 15 
IF (KTR.EQ.O) GO TO 13 
DO 12 KK*1,KTR 

12 IF (K.EQ.IABS( I SAME ( KK ) ) ) GO TO 15 

13 CONTINUE 

LET TT= COEFF EQ,X(K) 

1 S = 1 

LET TRUTH*MATCH 1D,SS,TT 
IF (TRUTH) GO TO 14 
LET TRUTH* MATCH ID,S$,-TT 
IF (.NOT. TRUTH) GO TO 15 
IS = -1 

14 KTR1=KTR1*1 
ISAME(KTR1)=K*IS 

15 CONTINUE 

16 CONTINUE 
MG=MG*1 

IF (MG.GT.MAXG) GO TO 20 
WRITE (6,24) 


D 56 
0 57 

D 58 
0 59 

0 60 
0 61 
0 62 
0 63 
D 64 
D 65 
0 66 
0 67 

0 68 

0 69 


0 72 

0 73 

D 74 
D 75 
D 76 
0 77 

0 78 

0 79 

D 83 
0 81 
0 82 
0 83 

0 84 

0 85 

D 86 
0 87 

D 88 
0 89 

D 90 
D 91 
D 92 
0 93 

D 94 
D 95 
D 96 
D 97 
0 98 

D 99 
D 100 
D 101 
D 102 
0 103 
D 104 
D 105 
D 106 
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D 107 


WRITE (6,211 MG 
PUNCH 21, MG 
CALL SOUT(SS,3, IH+) 

WRITE <6,231 MG, MG, MO 0 109 

PUNCH 23, MG, MG, MO 

LET DD=G(MG> D 119 

L1=KTR*1 0 112 

DO 17 J-L1,KTR1 D 113 

M^lABSt ISAMEI J>> D 114 

S=ISIGN(1,ISAHE<JI) 0 115 

LET EQ=EQ-S*SS*X(M> 0 116 

LET EQ» EXPAND EQ+$*D0*X1 M) D 117 

17 CONTINUE D 118 

ERASE DD,TT * SS D 119 

KTR=KTR1 0 120 

IF (KTR.EQ.NXS) GO TO 19 D 121 

18 CONTINUE 0 122 

19 CONTINUE D 123 

ERASE RR D 124 

RETURN D 125 

20 WRITE (6,26) D 126 

LET EO-EQSAVE 0 127 

RETURN D 128 

ENTRY COMOUT 0 129 

MO=MAXO( MO, 1 ) 

PUNCH 25, MG, MO D 139 

RETURN 0 131 

C 0 132 

C 0 133 

21 FORMAT ( 7X,2HG(,I4,2H>-) D 134 

22 FORMAT ( 7X,6HDEN0M< , 14, 2H) =) 

23 FORMAT ( 7X, 2HG( ,14, 4H)=6(,I4,8H)/DEN0M(,I4,1H1) 

24 FORMAT ( 1HX ) D 137 

25 FORMAT (6X, lOHCOMPLEX G(, 14, 8H>,OENOM( ,1 4, 1H) ) 0 138 

26 FORMAT (105HK**** WARNING THE NUMBER OF SUPER FUNCTIONS IG’SJ EXCE D 139 

1EDS THE MAXIMUM NUMBER OF FUNCTIONS (G'S) PERMI TTED , /, IOX, 43HTHERE 0 149 

2FORE SUBSTITUTION PROCESS STOPPED ***+> 0 141 

27 FORMAT (7X,2HG(, 14, 9H )=( DENQMl , I4,12H)-l.) /DENOMt ,I4,1H) ) 

END D 142- 
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Subroutine INSTRN 
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SUBROUTINE lNSTRN I FINISH) 

COMMON /STRING/ NUMBX ( 200 > . MM 
01 MENS I ON AS7RNG(14), INUM<4) 

logical entr 

DATA ENTR/, FALSE* / 

LOGICAL FINISH 

DATA BLANK. DOLRSN . RGHTPR . COMMA/ 1H t 1H$ , 1H) f IH. / 

DATA INUM/4+0/ 

DATA MBL/0006060606060/ 

FINISH*. FALSE* 

11*1 

MM=1 

IF (ENTR) GO TO 2 
ENTR** TRUE* 

1 READ <5, 10) { ASTRNGt 14) 

1 = 1 

2 CONTINUE 

STRNG=OR(FC( ASTRNG, D.MBU 
IF ( STRNG+ E Q • BL ANK ) GO TO 8 
C NOT A BLANK. SEARCH FOR A NUMBER 

NUMBER=LGR( 30. STRNG ) 

IF (NUMBER.LT*0.0R*NUMBER.GT*9) GO TO 3 
C FOUND ONE OF THE INTEGERS N OF THE X SUBSCRIPT X(NNNN) 

IF (IUGT.4) GO TO 9 
INUMC 1 1 »=NJM8£R 

n=n+i 

GO TO 8 

3 IF ( STRNG* NE*DOLRSN I GO TO 4 

C DONE WITH ALL INPUT STRINGS . TERM I NATE BUFFERS AND RETURN 

FI NI SH=*TRUE • 

GO TO 7 

4 IF ( STRNG ,NE* RGHTPR) GO TO 6 

C TERMINATE NUMBER BUFFER AND CONSTRUCT NNNN OF X(NNNN) 

11*11-1 
ISUM=0 
I TEN=1 
00 5 J*ltl! 

L*I I-J+l 

I SUM* 1SUM+1 NUM( L) * I TEN 
IT EN= ITEN+10 
INUM( L ) =0 

5 CONTINUE 

IF 1 1 SUM.GT *600 ) GO TO 9 
NUM8X ( MM ) =1 SUM 
MM=MM+1 
11=1 
GO TO 8 

6 IF ( STRNG .NE.COMMA) GO TO 8 

IF (PREV*N6. COMMA) GO TO 8 

C TERMINATE CURRENT STRING BUFFER 

7 CONTINUE 
MM=MM-i 
1*1 + 1 
RETURN 

8 PREV= STRNG 
1=1+1 

IF ( I *LE*8D ) GO TO 2 
GO TO i 


E l 
E 2 
E 3 
6 4 

E 5 
E 8 
E 6 
E 7 

e 9 
E ID 
E 11 
E 12 
E 13 
E 14 
E 15 
E 18 

E 17 
E 18 
E 19 
E 2D 
E 21 
E 22 
E 23 
E 24 
E 25 
E 26 
E 27 
E 28 
E 29 
E 3D 
E 31 
E 32 
E 33 

E 34 
E 35 
E 36 

E 37 
E 3B 

E 39 
E 4D 
E 41 
E 42 
E 43 
E 44 
E 45 
E 46 
E 47 
E 48 
E 49 
E 50 
E 51 
E 52 
E 53 
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o cr* 


WRITE ERROR MESSAGES 
CONTINUE 
WRITE 16,111 
STOP 
C 

10 FORMAT 1 13A6t A2 I 

11 FORMAT UOX ,45HSUBSCR1PT TOO LARGE FOR THIS VERSION 
END 


E 5* 
E 55 
E 55 
E 5T 
E 58 
E 60 
E 61 
E 62 
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Subroutine SOUT 
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ui^o oj ro ooooo 


^ U0KLJU f 1 nc. 5UU1 lAfKUUCt riMiir/ 

DIMENSION DUT( 12) 

SVMARG X 

DATA BLANK/ 1H / 

1 CODE-1 PAINT 
IC0DE=2 PUNCH 
I C0DE=3 PRINT AND PUNCH 


F 2 
F 3 

F 4 
F 5 
F 6 
F 7 
F 8 


C= PRNTP 
CONTINUE 

LET Q= BCDCON X,0UT t 12 
CALL REMOVE (OUT) 

IF i l CODE • E Q « 1 ) GO TO 2 

punch 4 , (ourn >» 1 = 2 , i2> 

IF < I CODEo E Q« 2 ) GO TO 3 
WRITE<6,5) C»C0UTtll«l»2tl2> 
C^BLANK 

IF (Q.NE.O. I GO TO I 
RETURN 

FORMAT I5X, 1H*, 11A6I 
FORMAT (Altl9X f llA6) 

END 


F 9 
F 10 
F 11 
F 12 
f 13 
F 14 


F 16 
F 17 
F 18 
F 19 

F 21- 
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Subroutine REMOVE 














SUBROUTINE REMOVE (X) 1 

INTEGER XU) 2 

LOGICAL NBR 3 

DATA AAST , AMINOS, AE , ALP, ARP, ADP, AIERQ/6H*0000D» 6H-00000 , 6HEOOOOO, S 4 

1 HI 00000 * 6H| 00000, 6H, 00000 , 6HOOOOOO/ , MSW/1 / 5 


DATA ADLR, K/ 6H$00000, 1 / 

LMAX=X( l)+6 
L“7 

1 GO TO (2,3, 4, 6), K 

2 IF <FC(X,L)*NE. AAST) GO TO 11 
GO TO 12 

3 IF (FC(X,LUNE. AAST) GO TO 11 
GO TO 13 

4 IF IFCIXfL) .EG. ALP) GO TO 14 

5 IF (NBR ( X, L ) ) GO TO T 
GO TO 11 

6 IF (FC(XfL) #NE* AM TNUS I GO TO 11 
L=L+i 


MS W=2 
GO TO 5 
T t=ua 

IF (NBR(X,LH GO TO 7 
IF (FC(X,L)*NE« ADP ) GO TO 11 
L*L + 1 

IF ( FC( X,L) *NE* AZERO) GO TO 11 
L-L + l 

GO TO (8,9) , MSW 

B IF (L.EG.LMAX+1 > GO TO 10 

IF ( NBR (X,L)»QR»FC(X,L)«EQ«AE) 
GO TO 10 

9 MSW*1 

IF ( L*EQ'LMAX+1 ) GO TO 10 

IF ( FC (X, L) • NE* ARP ) GO TO 11 

10 CALL DELETE (X,L-l) 

CALL DELETE (X,L-21 

11 K*1 

GO TO 15 

12 K=2 

GO TO 15 

13 K*3 

GO TO 15 

14 K-4 

15 L-L + l 

IF ( L« LE*LMAX ) GO TO 1 

IF ( FC ( XtLMAX ) • NE. ADLR ) RETURN 

K* 1 

CALL OELETE <X,LMAX> 

RETURN 

ENO 


6 

7 

8 
9 

10 
11 
12 

13 

14 

15 
IS 

17 

18 

19 

20 
21 
22 

23 

24 

25 

25 
27 

26 

GO TO 11 29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 
49- 


\ 

\ 

\ 

\ 

\ 
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FUNCTION FC IX, L) i 
DIMENSION X(l» ‘ 
DATA MASK/0770000000000/ 3 
LM1*L-1 t 
IW=LMl/6+l ' 
FC=ANDI ALSI 6*M0DI LM1, 6) ,X ( IW) ) ,MASK ) 6 
RETURN 7 
ENTRY FCltX, LI 8 
FC S ANDI ALSI L,X ) , MASK1 9 
RETURN 10 
END 


SUBROUTINE DELETE tX,L» 1 

DIMENSION XILI, ZER0I6I, BLANK! 61 2 

DATA ZERO/0007777777777, 0770077777777,0777700777777, 0777777007777, 3 

107777777700 77, 07777777777007, BLANK/0600000000000, 0006000000000,000 A 

20060000000 , 0000000600000 , 0000000006000 , 0000000000060 / 5 

LM1=L-1 6 

IW=LM1/ 6+1 7 

IP=MDDILM1,6H-1 8 

X ( 1 W) =OR( AND( X t IN I, ZERO! IP II, BLANK! IP I I 9 

RETURN 

END 11 


LOGICAL FUNCTIONNBRIX.L) 

REAL XIII 
N«LGRC3D,FCIX,L)t 
NBR=. FALSE. 

IF I 1N.GE.0I.AND.IN.LE.9)) NBR=.TRUE. 

RETURN 

END 


SUBROUTINE COSIXI 
ENTRY sinixi 
ENTRY TANHIXI 
ENTRY ATANIXI 
ENTRY ALOGIXI 
ENTRY DUMP! XI 
ENTRY P0UMP1X) 
ENTRY FMCOIFIX) 
ENTRY EXPOMP(X) 
ENTRY FMCOMP(X) 
STOP 
END 
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Evaluation Program - EVAL II 
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nn M 


REAL ITMP.K, LENGTH \ 

INTEGER PLTTYP z 

COMPLEX NO.OO.Q 3 

COMMON /ARRAY/ BA( 14) , NNUM.NN 4 

COMMON /CEVAL/ T(10 ) ,Z (10 ) ,WN( 10) t K (10 ) » NEXP.NS IN, NCOS , NS INH, NCOSH 5 

1 ?N SORT 6 

COMMON /OMEGA/ W,D0,N0<5) 7 

COMMON /NONLIN/ EXP A( 5 ) ,SINA(5 ) ,COSA< 5 ) ,EXPB(5 ) , SINBI 5 ) ,COSB( 5 ) , EX 8 

1PC (5) tSINC (5) tCOSC (S) .EXP0I5) »S IND( 5 ) ,COSD( 5 ) ,S INHA< 5 1 ,CO$HA( 5 ) , SQ 9 

IRTA ( 5) , SINH8(5 ) ,C0SHB(5 ) , $QRTB< 5 ) , S INHC <5 ) « COSHC< 5 ) * SQRTCI 5 ) . SINHD 10 

2(5),COSHD(5) 11 

COMMON /CLPLOT / XPEN, YPEN ,NX ,NY , IPEN.XLABEL ( 10) * YLABEL ( 10 ) 12 

DIMENSION Yl(604,5), Y2(604»5)» KKK(6), PC14)# X(604> 13 

DIMENSION NDENI5) 14 

DATA LENGTH/O* / I 5 

DATA L0GPL0/6HL0GPL0/ l 6 

DATA CONVT /5T. 2957795/ I 7 

DATA YDB/3HDB./ 1® 

DATA (KKK ( J) ,J*2»6) ,P<1) , (P(L ) ,L*3, 14 ) /50, 1 , 1,0, 601, 3. , 2*0. , 10., 2* 19 

10. ,10.. 5*0. ,90./ 20 

DATA TWOP I / 6. 2 8318 53/ 21 

OATA XLAB/lHF/ 22 

DATA YLAB/3HABS/ 23 

DATA ZLAB/5H ANGLE/ 24 

DATA NK,NWN,NT,N2/4*0/ 25 

NAMELIST /INPUT/ NEXP,NSlN,NCOS t EXPA,SlNA,COSA,EXPB,SlNB»COSB,EXPC 26 

1, SINC.COSC, EXPO, SINO, COSO, NSINH, NCOSH, NSQRT,SINMA,CO$HA,SQRTA,SINH 27 

2B .COSHB , 5QRTB, SINHC .COSHC »SQRTC,$ INHO.COSHD »NK*K,NT ,T » NZ» Z » NWN, WN» 28 

2FSTART, FEND, DELTAF, LENGTH 29 

C READ THE NUMERATOR AND DENOMINATOR INFORMATION SUPPLIED BY REOUCE II 

RE AO (5.12) NNUM, IK, (NOEN (I ), 1*1, IK) 30 

C READ THE NAMELIST INPUT. NOTE THAT THE PROGRAM ALWAYS RETURNS HERE 

1 READ 15, INPUT) 31 

WRITE (6,18) 32 

IUPPER-1 33 

C CHECK FOR A PLOTTING REQUEST. 

IF (LENGTH.EQ.O.) GC TO 2 34 

IUPPER-2 35 

C PLOTTING TO BE DONE. REAO IN PLOT TITLE ANO PLOT TYPE 

READ (5,30) BA 36 

WRITE (6,31) BA 37 

READ (5,30) PLTTYP 36 

CONTINUE 39 

EXAMINE THE VARIOUS NAMELIST VARIABLE TO GET THE INPUT IDENTIFIED , 

ON THE PROGRAM OUTPUT LISTING ' 

IF (NEXP.NE.O) WRITE (6,19) (l.EXPA(I), EXPB( I) ,EXPC( I),EXPD(I),I*1 40 

l.NEXP) *1 

IF ( NSIN.NE.O) WRITE (6,20) ( I ,SINA( I ) ,SINB ( I) ,S INCI I ) , SIND(I) , 1*1 42 

I. NSIN) A3 

IF (NCOS.NE.O) WRITE (6.211 ( I ,COSA( I ) ,C0S8( I) ,COSC( I ) ,COSO( I ) » 1*1 44 

1 .NCOS ) 45 

IF (NSINH.NE.O) WRITE (6,221 ( I »S INHA1 I ) ,S INH8( I ) , SINHC ( I ) , S1NHD( I 46 

l),I-l,NSINH) 47 

IF ( NCOSH. NE.O) WRITE (6,23) ( I «COSHA( I ),COSHB( I ) , COSHC (1 ),COSHD( I 48 

II, 1*1, NCOSH) 49 

IF ( NSQRT.NE.O I WRITE (6,24) ( I, SORT A ( I ),SQRTB< I), SQRTC(I ), I»1,NSQ 50 

1RT) 51 

IF (NK.NE.Q) WRITE (6,25) ( I ,K( I ) , 1*1, NX) 52 

IF (NT.NE.O) WRITE (6,26) < I ,T( I ) , 1*1, NT) 53 
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IF (NZ.NE.C) WRITE <6« 27) . f I t Zfl) » 1-1 »NZ) 




54 


IF (NWN.NE.O) WRITE (6*28) C I t WN( I ) » 1*1 *NWN ) 




55 


WRITE (6«13) UfNNUFtNOENd), 1*1, IK) 




56 


WRITE (6*29) F START* FEND* CELT AF 




57 


WRITE (6*14) 




58 


IF (DELTAF* EQ.O.O) GO TO 1 




59 

c 

START THE MAIN PROGRAM LOOP* IUPPER *1 IF THERE ARE 

NO 

PLOTS 



c 

2 IF PLOTS ARE 

TO 

BE MADE 




DO 8 NTYPE*1 *1 UPPER 




60 


IF (NTYPE.EQ.2) GO TO 3 




61 


DF-OELTAF 




62 


GO TO 4 




63 

3 

OF* ( PEND-F START ) /600. 




64 

A 

J*0 




65 


F*FSTART 




66 

5 

W*TWOPI*F 




67 


IF (J.GE.60A) GO TO 8 




68 
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C 

THE CALL TO COEFF EVALUATES THE TRANSFER FUNCTION AT 

S*I*W 




CALL COEFF 




70 


Xf J)*F 




71 


00 7 1*1,1* 




72 


G*N0U)/00 




73 


RTMPasRE At ( Q) 




74 


I TMP^AI MAG(G) 




75 


Yiiu,n*o. 




76 


Y2<J*I>*0. 




77 


ZCB*0. 




78 


IF < RTMP.EQ.Q. » AND* ITMp.EQ.O* 1 GO TO 6 




79 


Y I ( J * I ) *C ABSI Q) 




80 


ZCB*20**ALOGIOIY1( J*I) ) 




81 


Y 2 ( J * I) *A TAN2 (1 TMP * RTMP > *CONVT 




82 

6 

CONTINUE 




83 


IF (PLTTYP*EQ. LOGPLG. AND* NTYPE.EQ.2) Yl(J t I>*ZDB 




84 


IF (NTYPE.EC.2) GO TO 7 




85 


IF ( I .EQ.1 • AND. 1K.GT.1) WRITE (6*15) 




06 


WRITE ( 6*16) I*F,RTMP,ITMP,Yl(J f I),Y2( J,!)* ZDB 




87 

1 

CONTINUE 




88 


F*F+DF 




89 


IF (F.LE.FENO) GO TO 5 # 




90 

8 

CONTINUE 




91 

C 

CONTINUE ON WITH THE PLOTTING IF IT IS CALLED FOR 






IF (LENGTH. EQ.O.O) GO TO l 




92 


K KK ( 6 ) * J 




93 


PC2UAINT(LENGTH+ # 5> 




94 


NX*-1 




95 


XL ABEL ( 1 )»XL AB 




96 


OC ii I *1 * IK 




97 


P(6)*0. 




98 


P(?MO. 




99 


NY *3 




100 


IF (PLTTYP.EQ.LOGPLG) GO TO 9 




101 


KKKU)«4 




102 


YL ABEL ( 1 )*YLAB 




103 


GO TO 10 




104 

9 

CONTINUE 




105 


KK K 1 1 ) »2 




106 


IF (I.GT.U KKK(l)*-2 




107 


YLABEL ( 1 )*YDB 




109 

1C 

CONTINUE 




lG r 
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n n 


NA-NDENII) 

C MAKE THE AMPLITUDE RATIO PLOT 

CALI LOGPLT IX ,Y1 < 1 ,1) ,KKK,P> 
WRITE (6,17) 

IF IPLTTYP.EQ.LOGPLOl KKK<11«-2 

PI6)*-200.00 

P(7)*2GQ.OO 

YLABELIli-ZLAB 

NY*5 

C MAKE THE PHASE ANGLE PLOT 

CALL LOGPLT IX,Y2U »1 1 ,KKK,P> 
WRITE (6,17) 

II CONTINUE 

GO TO l 


12 

12 

1* 


15 

1C 

n 

ie 

19 

20 
21 
22 
23 
2* 

25 

26 
27 
2£ 

29 

30 

31 


TRANSFER FUNCTION* </»5ZX, 12* 9X*2HXI *i!2* AH 


FORMAT ( 15 ) 

FORMAT I 50X* 27HNUM8ER 

^ FORMA 7^ I //»3X, 6HNUMBER* 5X »9HFREQUENCY * 13X*9HREAL PART»9X, 14HIKAGIN 
1 ARY PART.6X* l*H ABSOLUTE VALUE *6X, 11HPHASE ANGLE*9X» 16HAMPLITIDE IN 
2 OB.) 

FORMAT I1HK) 

FORMAT <I5*6F20.6) 

FORMAT (10X.9HPL0T MADE) 

FORMAT ( 1H1 * ft *36H FORTRAN EVALUATION OF FORMAC OUTPUT) 

FORMAT (1H0»3X»1HI *5X*5HEXPA ,10X,5HEXPB ,10X,5HEXPC »10X,5HEXPD * 

^F0RMAT^^IH0»3X »1HI * 5X* 5HSINA ,10X,5HSINB ,10X,5HSINC »10X,5HSIN0 * 

Format* ciho*3x»ihi »sx»5hcosa ,iox*5hcosb ,iox,5hcosc *iox,5hcosd ♦ 

Format 4 C1Ho!3X» 1HI*5X»5HS INHA, 10X,5HSINHB,10X,SHSINMC, 10X* 5HSINH0* 

1 FORMAT^ 1H0,3X #1HI , 5X, 5HC0SHA * 10X »5HC0SHB, 10X, 5HC0SHC * 10X , 5HCOSHO, 

^FORMAT' 4 '( 1H0^3X»1HI * 5X»5HS0RTA»i0X»5HSQRTB»10X,5HSQRTC»/ »II5,3E15.5 
111 

FORMAT llH0,3X,lHi,5X,lHK»/,< I5,E15.5) ) 

FORMAT <1H0,3X,1HI »5X,1HT,/, I I5»E15.5) ) 

FORMAT (1H0,3X,1HI »5X,1HZ,/,I 15,615.5) ) 

FORMAT (1H0,3X,1HI ,5X,2HMN,/ , 115, E15.5 )1 

FORMAT |//,8H FSTART-,F15.5,5X,5HFEN0«»F15.5,5X,7H0ELTAF*,F15.5) 
FORMAT ( 13A6 ,A2 I 

FORMAT (3AHOCALCOMP PLOTS HAVE BEEN REQUESTED, //,5X,13A6, A2) 

END 


110 

111 

112 

113 

11 * 

115 

116 

117 

118 

119 

120 
121 
122 
123 
12 * 

125 

126 
127 
12B 

129 

130 

131 

132 

133 
13* 

135 

136 

137 

138 

139 
1*0 
1*1 
1*2 
1*3 
1 ** 
1*5 
1*6 
1*7 
1*B 
1*9 

150 

151 

152 
153- 
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Subroutine COEFF 


SUBROUTINE COEFF 

C PLT REDUCE II OUTPUT IN THIS SUBROUTINE 

REAL K 

£****« ******************************************************* *********** 
C* * THE COMPLEX STATEMENT PUNCHED BY REDUCE II FOLLOWS THIS CARD * * * 
C *********************************************************** ************ 

LOGICAL OQNE 
DATA DONE/. FALSE./ 

COMPLEX NOtOOt E XPFX t SI NFX yCOSFX «S INHFX 9 COSHFXf SQRTFX* $ 

COMPLEX EXPF(5> ,$INF<5) ,CO$Fm,SINHF<5 >,C0SHF(5)t $QRTF(5> 

COMMON /CEVAL/ T < 10 > * Z < 10 > « WN UO ) • K( 10 I t NEXP,NS IN, NCOS, NS INH f NCOSH 
IfNSQRT 

COMMON /OMEGA/ W*00»N015> 

IF i NEXP.EQ* 0) GO TO 2 
DO 1 I»X*N EXP 

1 EXPFm-fEXPFX( I) 

2 IF (NSIN.EQ.O) GO TO 4 
00 3 1*1 »NSIN 

3 S I Nf t I ) * $ I NF X t I) 

A IF { NCOS. EQ.O) GO TO 6 

DC 5 I =* I, NCOS 
5 COSF < I ) »C0 SFX{ I I 

£ IF (NSINH.EG.Of GO TQ 8 

DO T 1*1 1 NS INH 

7 S INHF tl )*S INHFX { I ) 

8 IF ( NCOSH. EC.O ) GO TO 10 

00 9 f ®1 * NCOSH 

S CCSHFU MCOSHFxm 

IQ IF INSORT. EG. 01 GO TO 12 

DC 11 1*1 t NSQRT 

11 5QRTFU WSQRTFXU) 

12 CONTINUE 
S*CMPLX (0.0 *W) 

C****<******************************************************^*********** 

c * * * THE USERS • S FUNCTIONS (G(Sl i S) ARE INSERTED AFTER THESE CARDS * * 
C* * * THE G f S MUST BE FOLLOWED BY THE SUPER G«S,THE NUMERATOR f AND * * 
C* * * DENOMINATOR CARDS REDUCE II PUNCHED * * * 

£****« ************** ******************************** ******************** 
RETURN 
END 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 
36 
39 
40 ~ 
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Function EXPFX 


COMPLEX FUNCTION EXPFXCI) l 

COMPLEX X»Y 2 

COMMON /OMEGA/ W,D0,N0<5) 3 

COMMON /NONLIN/ EXP AC5 ) ,$ INAC5 1 »COS AC 5 1 ,EXPB 1 51 »SINB( 5 ) .COSBC 5 1 ,EX 4 

lPC(51fSINC(51*COSC(51fEXPD(5l*SINQ(51« COSDC 5 I »S INHAC 5 1 *COSHA{ 5 1 * SQ 5 

<RTAC5),SINHB(5) »COSHB C5 1 , SQRTB 1 5 ) ,S INHC C5 ) ♦ COSHC < 5 1 , SQRTCC 5 > , S INHD 6 

315) «COSHD {51 7 

X*CMPLX (EXPC ( I >-EXPA(I)*W**2,EXPBC I1*W) 8 

Y«EXPDCI)*C$QRTCX) 9 

EXPFX=CEXP( Yl 10 

GC TO 1 11 

ENTRY SINFXU) 12 

X- CMPLX(SINCm-SINAm*M**2,SIN8f I1*W) 13 

Y*SINDCI1*CSQRTCX) 14 

EXPFX-CSIN(Y) 15 

GO TO 1 16 

ENTRY CCSFXCI1 17 

X*CMPLX (COSC ( 1 1 -COSAC l 1*W**2,C0SBC I >*VO 18 

Y»COSDll)*CSQRTCX) 19 

EXPFX=CCOS ( Y ) 20 

GC TO 1 21 

ENTRY SINHFXCIl 22 

X*CWPLXCSI NHC Cl 1— SINHACI)4W*+2tS I NH8 (!) *W ) 23 

Y “ S tNHD C I) *C SQR T C X ) 24 

X1«REAL (Y ) 25 

Y1=AIMAGCY) 26 

EXPFX=CMPLXCSINHCX1)*C0SCY1)»C0SHCX1)*SINCY1)1 2? 

GC TO 1 28 

ENTRY COSHFXU) 29 

X a CMPLX CCOSHCC I ) -COSHA (I) *W**2»C0SHB < 1 1 +W 1 30 

Y*COSHD 1 1 )*CSQRT C X) 31 

XI- REALCY) 32 

Y1*AIMAG( Y) 33 

EXPFX=»CMPLXCCOSHCX1)*COSIY1) .S I NH C XI ) *S INC Y 1 }) 34 

GC TO 1 35 

ENTRY SQRTFXCI) 36 

X*CMPLX CSQRTCCI1-SQRTACI1 *W**2 » SQRTBC I)*W) 37 

E XPF X=C SQRT C X 1 36 

l RETURN 39 

ENO 40- 
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Subroutine CALTIT 


SUBROUTINE CALTIT 

COMMON /ARRAY/ BAtlA) ,NDEN,NNUM 

DIMENSION FUNC(A) 

DATA FUNCU W20HTRANSFER FUNCTION X / 

ACsNOEN 

AN*NNUM 

CALL SYMBOL 1-1.5 , l. ,. 15, BA, 90. »80 I 
CALL SYMBOL (-1 .0, 1. , . 15, FUNC ,90. ,19 > 
CALL NUMBER (-. B5 , 3.*,. 100, AD ,90. ,-l ) 


CALL SYMBOL (-1. ,3.52, .15 ,ZH/X,90. , 2 > 10 

CALL NUMBER <-. 85, 3. 70, .100, AN, 90. ,-l ) 11 

RETURN 12 

END 13- 


Block Data 


BLOCK DATA 

COMMON /CLPLOT/ XPEN, YPEN ,NX , NY, I PEN, XL { 10 » , YL< 10 1 
DATA XPEN , YPEN , I PEN/2*0. ,-3/ 

END 
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Subroutine LOGPLT 


StSROUT IN 6 LOGPLT < XDOWN, YACROS ,KKK f P ) 1 

2 

THIS LOG PLOTTING PROGRAM WAS WRITTEN BY MR. BERT HENRY OF 3 

LEWIS RESEARCH CENTER. IT WAS MODIFIED BY «R. JERRY LENHART AND 4 

MR. JOHN RIEHL ALSO OF LEWtS RESEARCH CENTER 5 

6 

THIS IS A GENERAL ROUTINE FOR MAKING LOG-LOG ANO SEMI-LOG CALCOMP 7 

PLOTS ON PLAIN PAPER. 8 

9 

this subroutine utilizes the same conventions as calplt except for io 

KKKII1 ,pm ,P<4),PI6> , AND PIT) AS INDICATED BELOW. 11 

THE ROUTINE EXAMINES KKK(l) TO DETERMINE WHICH AXIS WILL 12 

REQUIRE A LCG SCALE t THEN EXAMINES THE APPROPRIATE DATA IB 

TO DETERMINE THE NUMBER OF CYCLES REQUIRED. IF THE NUMBER 14 

OF CYCLES DOES NOT EXCEEO 7 OR 15 < Y OR X AXIS RE- 15 

SPECTIVELY » , THE GRIDS ARE DRAWN 16 

PERPENDICULAR TO THE LOG AXIS AT EVEN VALUES IN EACH CYCLE 17 

IF MORE THAN 7 OR 15 CYCLES ARE REQUIRED* AN APPROPRIATE 18 

MESSAGE IS PRINTED AND CONTROL IS RETURNED TO THE CALL- 19 

ING PROGRAM. 20 

21 

FOR THE LINEAR SCALE. THE ROUTINE WILL COMPUTE THE SCALE IF ANO ONLY 22 

IF BOTH PIBI AND PC41, OR P( 6 ) ANO P(7», ARE SET TO 0.0 23 

FOR THE X OR Y AXIS RESPECTIVELY. GRID LINES WILL BE 24 

DRAWN PERPENOICULAR TO THE AXIS AT THE TICK MARKS WHICH 25 

ARE CONTROLLED BY PIBI. 26 

27 

A CESCRIPTION CF THE SUBROUTINES ARGUMENTS FOLLOWS-* 2B 

XDCWN IS THE NAME OF THE ARRAY CONTAINING THE VARIABLE TO BE PLOTTER 29 

AS THE ABSCISSA 30 

YACRCS « THE NAME OF THE ARRAY CONTAINING THE VARIABLE TO BE PLOTTED 31 

THE ORDINATE 32 

33 

KKK ( l ) = 1* IF Y AXIS ONLY IS LOGARITHMIC 34 

KKM1) 4 2 , IF X AXIS ONLY IS LOGARITHMIC 35 

KKK ( 1 1 - 3, IF BOTH X ANO Y AXES ARE LOGARITHMIC 36 

KKKIU - 4, IF BOTH X AND Y AXES ARE LINEAR 37 

*«♦* IF KKKU> IS NEGATIVE THE APPROPRIATE NON-LINEAR DATA IS 38 

ASSUMED TO BE IN LOG BASE 10 FORM. IF POSITIVE THEN THE DATA 39 

WILL BE TRANSFORMED TO LOG BASE 10 AND LEFT IN THIS FORM. 40 

KKK( 2 ) * USED TO CONTROL SYMBOLS PLACED ON DATA POINTS 41 

KKM 2 I 0* FOR ONLY A LINE PLOT 42 

KKK C 2 1 1.2.3, ETC, WILL PRODUCE A SYMBOL AT EVERY DATA POINT, OR 43 

EVERY KKKI2I DATA POINTS 44 

A NEGATIVE KKK(2 ) WILL SUPPRESS THE LINE BETWEEN DATA POINTS 45 

KKK ( 3 1 * THE NUMBER OF CURVES TO BE PLOTTED IKN 1 46 

KKM4) = THE REPEAT CYCLE IN CASE VALUES ARE IN A MIXED ARRAY 47 

KKK ( 5 > = THE DESIRED SYMBOL 4 B 

KKKI 6 I - KKKIKN+5I « THE NUMBER OF POINTS IN EACH CURVE 49 

P < 1 » * 1.0, OPTION FOR OUPX, AS USED IN PLOTMY 50 

P< 1 ) * 2.0, OPTION FOR OUPY, AS USEO IN PLOTMY 51 

P(l) = 3.0, OPTION FDR NODUP, AS USED IN PLOTMY 52 

P ( 2 ) * THE MAXIMUM WIDTH OF THE PLCTTEO LINE 53 

*** **NOTE***** 54 

IF ABSCISSA VALUES ARE LINEAR, THEN 55 

P(2) * MINIMUM VALUE OF ABSCISSA ARRAY OF DATA 56 

P(<) * MAXIMUM VALUE OF ABSCISSA ARRAY OF DATA 6 ? 

IF ABSIKKKUN * 1 AND EITHER P<31 OR P(4> IS NON-ZERO, THEN THF X 5 R 
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o r» n n o o 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SCALE WILL NOT BE COMPUTED BY THE PLOT ROUTINE* 

P(M * THE MAXIMUM HEIGHT OF THE PLOTTED 
MUST NOT EXCEED 10. INCHES. 

♦ M**NOTE***** 

IF ORDINATE VALUES ARE LI NEAR f THEN 

PI*) * THE MINIMUM VALUE OF THE ORDINATE ARRAY OF DATA 
pn) n the MAXIMUM VALUE OF the ORDINATE ARRAY OF DATA 

IF ABS (KKK ( 1 ) ) * 2 AND EITHER P<6) OR P 1 7 1 IS NON-ZERO, THEN THE Y 
SCALE WILL NOT BE COMPUTED BY THE PLOT ROUTINE. 

♦♦***ncte***** 

pm^DIv! °5lLt A CAUSE N A A TICK M MARK TO BE PLACED AT EVERY 10.0/DIV IN 

the coordinates of the starting point of the abscissa 

pill l - P( 12 ) »* THE COORDINATES OF THE STARTING POINT OF THE ORDINATE 

Pf 13) * THE ANGLE IN (DEGREES! OF THE ABCISSA AXIS 

PU4) * THE ANGLE IN (DEGREES! OF THE ORDINATE AXIS 

PC15! =* NON-ZERO IF DATA IS TO BE ROTATED PER P(13) AND P(14! 

NOTE******** 

THE FOLLOWING CARO IS NEEDED IN THE CALLING PROGRAM 

COMMON /CLPLOT/XPEN,YPEN,NX f NY,IPEN,XLABEL(IOJfYLABEL( 10! 
uy to c«* 

XPEN AND YPEN GENERATES THE INCREMENTS NECESSARY TO HOVE FROM CUR 

POSITION TO XPEN, YPEN NORMALtYaO.O 

IPEN*3, FOR PEN UP 

lPEN--a’oR-2 ^PROVIDES A NEW REFERENCE POINT AT (XPEN, YPEN! RELATIVE 
THE CURRENT REFERENCE POINT 

XLABEL * THE NAME OF THE ALPHABETIC ARRAY OF DATA FOR THE ORDINATE AX 
YL ABEL * THE NAME OF THE ALPHABETIC ARRAY OF DATA FOR THE ABCISSA AXI 
NX* THE NUMBER OF CHARACTERS IN THE ORDINATE AXIS TITLE, A MINUS NX P 
THE ANNOTATION ON THE CLOCKWISE SIDE OF THE AXIS INSTEAD OF THE 
COUNTER CLOCKWISE SIDE 

NY* THE NUMBER OF CHARACTERS IN THE ABCISSA AXIS TITLE, A MINUS NY P 
THE ANNOTATION ON THE CLOCKWISE SIDE OF THE AXIS INSTEAD OF THE 
COUNTER CLOCKWISE SIDE 

COMMON /CLPLQT / XPEN , YPEN ,NX,NY, I PEN, XLABEL ( 10 ) , YLABEL (10 ) 

COMMON /SPECL/ TEST ,ORG$ET,$PASET 
COMMON /XCPIDX/ NBLCNO* IRUNO»CP 10(81 
DIMENSION SAVARA14! 

DIMENSION XOOWN(l), YACROSIll, KKK ( 1 ) , P(l) 

EQUIVALENCE (STR,ISTR! 

DATA RUNHES/6HRUN NC/ 

OATA BLKMES/6HBLK NO/ 

DATA INITAL/50/ 

OATA STR/OSCOOQOOOOCOO/ 

OATA SPACER/8./ 

IF ( INI TAL.EQ.50! GO TO 3 
IF (P(11.EQ.4.0! GO TO 2 
IF ( IPEN.NE.(-3)1 GO TO 5 
IF (XPEN.EQ.STR1 XPEN=0. 

IF (YPEN.EQ.STR1 YPEN*0. 

GO TO 4 

CALL PLOT ( XPEN, YPEN, I PEN! 

RETURN 
XTPEN-0.0 
YTPEN»0.0 
I TPEN=-3 
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XSPAC=SPACER 119 

IF IPm.6Q.STR) P( 1)*3.0 120 

IF <P(2).EQ.STR) P< 2 ) * 10.0 121 

IF <P(3).EQ.STR) Pt 3)=0. 122 

IF IPI4J.EQ.STR) PI4)*0. 123 

IF I P ( 5 ) . EQ. STR ) PI 5)=10.0 124 

IF I P 1 6 > »EQ. STR ) PI 6)*0. 125 

IF I P (7 ) «EQ. STR ) P I 7 ) =0 . 126 

IF IPt8J.EQ.STR) P I 8 ) * 10 • 0 127 

IF IP(R).EQ.STR) P(9)*0. 128 

IF IP(IC).EQ.STR) PI10M0. 129 

IF (P(ll).EC.STR) PI1U-0. 130 

IF I PI 12 ) .EC.STR) P (12 )*0. 131 

IF (P(13).eC.STR) P I l 3 ) =0 • 132 

IF I P 1 14 ) .EC.STR) P (14 )*90.0 133 

IF 1 PI15 ) .EC.STR) P (15 )*0* 134 

IF (WKK(l).EQ. I STR) KKK(1)*1 135 

IF IKKK(2).EQ.I$TR) KKK 12 )*0 136 

IF (KKKI3J.EQ.ISTR) KKK(3)=1 137 

IF (KKK (4 ) .EQ. I STR ) KKK(4)*1 130 

IF ( KKK ( 5 ) . EQ. I STR ) KKK(5)*1 139 

IF (KKK(6).EQ. ISTR) KKK(6)*1 160 

IN ITAL-O 141 

CALL PLOT <XTPEN,YTPEN,ITP£N) 142 

FBLONO-NBLONO 143 

CALL SYMBOL <0.0, 1.0,0. 2 ,CP 10 1 1 1 ,90.0, 24 ) 144 

CALL SYMBOL <0. 5 ,1.0, 0.2, RUNMES ,90.0,6 ) 145 

CALL SYMBOL (0. 5 ,2. 5,0.2 , IRUN0,90.0,6 ) 146 

CALL SYMBOL (0.5,4.0,0.20,BLKMES,90.0,6 ) 147 

CALL NUMBER (0.5,6.0,0.20,FBLONC,90.0,-1) 148 

XTPEN*1C.0 149 

CALL PLOT (XTPEN,YTPEN,ITPEN> 150 

GC TO 1 151 

4 CALL PLOT ( XPEN ,YPEN, I PEN) 152 

5 NX1*-1 153 

FBLONO»NBLONO 154 

CALL SYMBOL I-3.5,0.0,0.15,BLKMES,0.0,6) 155 

CALL NUMBER 1-2. 5,0.0, 0.15,FBL0N0,0. 0,-1) 156 

YTPEN*0.0 157 

XTPEN-2.0 158 

ITPEN—3 159 

CALL PLOT <XTPEN,YTPEN,ITP£N) 160 

C CHECK FOR SCALING IN DUPX OPTION 161 

JX*1 162 

J Y*1 163 

IF (PO).NE.O.) GO TO 6 164 

IF (P(4).EQ.O.) JX»2 165 

6 IF IPI6J.NE.0.) GO TO 7 166 

IF (pm.EQ.O.) JY-2 167 

7 J* I ABSI KKK 111) 168 

I POPT*P 1 1 ) 169 

KN*KKK( 3) 170 

NPTS-0 171 

OC 8 1*1, KN 172 

8 NPT$*NPTS+KKK( 1*5) 173 

C TRANSLATE AND ROTATE OATA IF OESIREO 174 

IF IP(3).NE.3.) GO TO 9 175 

C IF IPI15I.EC.0. ) GO TO 63 176 

CALL TRNSFM I XDOWN, YACROS,NPTS,P( 13 ) , P ( 14)) 177 

C SETUP NUMBER OP POINTS FOR DUPX, AND DUPY 178 
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9 NSPTS-KN*KKK<6) 

N P T *KKK (61 

IF OCKK«6),EQ.O> GO TO 29 
GC TO (10,11,121,1 PCPT 
C SETLP FOR AXIS SUBROUTINE 
C CHECK FOR DUPX,DUPY OR NOOUP OPTION 

1C NP 1 -NPT 

NP2*NSPTS 
GC TO 13 

11 NPl^NSPTS 
NP2=NPT 
GC TO 13 

12 NP1-NPTS 
NP2=NPT S 

13 GC TO 114,19,18,14) ,J 

14 GC TO <15,161,JX 

15 XP I N^P ( 3 l 
DELX-?P?4I-P<3) >/P<2) 

GC TO 1? 

14 CALL SCALE < XDO WN , P <2 ) # NP 1 , KKK (4 J , P ( 0 ) , XMI N , DELX ) 

17 CALL AXGRID < P (9) f P (10 ) tXLABEl , NX , P(2 ), Pf 13 > tXMIN, DELX , P< 0 ) r P< 5 ) f P 
14141) 

IF (J.EG.4) GO TO 19 

18 CALL GRID (P‘il!!»PC12) , YL ABEL »NY , YACROS ,NP2 ♦ KKK <1) , P? 5 ) , P<14> ,P< 2 ) 
1,P<13>,YHIN,DELY,7> 

GC TO <24,24,23,24) ,J 

19 GC TO 120,21), JY 

2C YP IN-P ( 6 ) 

DUY-(Pm-PI6> )/P<5> 

GC TO 22 

21 CALL SCALE { YACROS , P ? 5 ) , NP2 , KKK (4 1* P i 8 > , YMI N f DELY ) 

22 CALL AXGRID < P (11 1 , P< 1 2 ) » YLABEL , NY , P ? 5 ) , P 114 1 , YMIN , DEL V ,P ( B > *P < 2 ) , 
1PC13) ,15) 

IF U.EQ*4> GO TO 24 

23 CALL GRID < P < 9 ) , P < 10 1 , XLABEL ,NX , XDOMN , NP1 , KKK? 1 > ,P ? 2 ) , P <1 3 ) , P ( 5 ) , P 
1 ( 14) ,XMIN,DELX) 

C SETLP FOR ENTERING LINE AND AXIS SUBROUTINES 

24 CONTINUE 
K2*0 

«J J J2-0 

DO 28 I ND“1 t KN 
K 3 ^ I ND 
K4-K2 

IKSYMB=KKK(5)+IND-1 
I S YP^2 

IF ? KKK ? 2 ) # EQ* 999} l$YM*KN+K3+5 
C SETLP FOR CHECKING OUPX,OUPY,QR NODUP OPTION 
C 

GC TO <25,26,27) ,IPCP7 
C DUP > CPTION FOR LINE 
C 

25 IF { KKK ( 6 ) « E Q* 0 ) GO TO 28 
K2 = K 2+KKK ( 6 ) 

CALL LINE ? XOOVIN, YACROS (K4+1 ) ,KKK < 6 > , KKK? 4 ) ,KKK ( ISYM) , IKSYMB, XM1N , 

ideix t ymin,oely» 

GC TO 28 

C DUP > OPTION FOR LINE 
C 

26 IF { KKK l 6 ) # E Q, 0 ) GO TO 28 
K2-K2+KKK?6) 
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CALL LINE (XD0WN(K4n> , YACROS ,KKM6 > ,KKK|4) .KKKUSYM), IKSYM8..XMIN, 239 

1DELX,YM!N,DELV) 240 

GC TQ 28 2^1 

C NOOLP OPTION FOR LINE 242 

C 243 

21 IF <KKK(K3+5l.EQ.O) GO TO 28 244 

K2=K2+KKK(K3+5) 245 

K1*K2 246 

JJ42-K2+2 247 

CALL LINE IXD0WN(K4+l),YACR0S(K4+l),KKK{K3+5),KKK<4)fKKKIISYM>,IKS 248 
1YP0,XMIN,OELX»YN1N,OELY> 249 

28 CONTINUE 250 

CALL CALTIT 251 

XPEN1=P (2 l+XSP AC 252 

IPEN=XPEN1 253 

XPEN1=IPEN 254 

YPEN=0.0 255 

I PEN=-3 256 

CALL PLOT (XPENlfYPEN»IPEN) 257 

XSPAC=SPACER 258 

29 RETURN 

END 259- 
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Subroutine GRID 


c tms subroutine will craw axis with log scale and its grid lines l 

SUBROUTINE GRID I X , Y ,BC0,NN, XX, N, M,S IZE , THET A, HGT.THT Y, XXM IN^XXDEL 2 

l.KCYMAX) 3 

c WHERE- X,Y IS THE COORDINATE Of THE STARTING POINT OF THE AXIS. A 

C BCD IS THE LOCATION OF ALPHA INFORMATION FOR AXIS TITLE. 5 

C NN IS THE NUMBER OF CHARACTERS IN ECO TITLE. A NEGATIVE 6 

C NN PLACES THE TITLE ON THE CLOCKWISE SlOE OF THE AXIS 7 

C LINE ANO VICE-VERSA B 

C XX IS THE NAME OF THE ARRAY IN THE AXIS DIRECTION 9 

C N IS THE NUMBER OF POINTS GIVEN IN THE XX ARRAY 10 

C M IS POSITIVE IF THE XX ARRAY IS TO BE CONVERTED TO LOG 11 

C BASE 10. FOR NEGATIVE M THE XX ARRAY VALUES WILL BE 12 

C TAKEN AS LOG BASE 10 VALUES 13 

C SIZE IS THE LENGTH OF THE AXIS TO BE PLOTTED IN INCHES 14 

C THETA IS THE ANGLE IN DEGREES OF THE AXIS MEASURED IN 15 

C COUNTER-CLOCKWISE DIRECTION FROM THE X AXIS. 16 

C HGT IS THE LENGTH OF THE OTHER AXIS IN INCHES 17 

C THTY IS THE ANGLE IN DEGREES OF THE OTHER AXIS MEASURED 18 

C COUNTER-CLOCKWISE FROM THE X AXIS. 19 

C 20 

DIMENSION XXIII, TI C ( 9 ) » XSPC<9), YSPCI9), VAL ( 5 ) 21 

DATA (TIC II ), I *1,9)/. 0457575,. 05 115 25,. 057992,. 0669467,. 0791813,. 0 22 

19691,. 1249387,. 1760913,. 30103/ 23 

DATA IVALUI ,I = 1,5)/B. ,6. ,4.,2.,1./ 24 

DATA RAD/. 0174532925/ 25 

K A T£ * l 26 

XMIN=l.E28 27 

XMAX*— XMIN 28 

DO 2 1*1, N 29 

IF (M.LT.OI GO TO l 30 

XX< I)=ALOGlC<XX<n ) 31 

1 IF < XX( D.GT.XMAX) XMAX=XX<I) 32 

IF < XX< D.LT.XMIN) XKIN=XXm 33 

2 CONTINUE 34 

KXKAX=XMAX 35 

XHAXK=KXMAX 36 

IF (XMAX.GT.XMAXK) XMAXK=XMAXK+1. 37 

KXM I N=XKI N 38 

XMINK=KXMiN 39 

IF (XMIN.LT.XMINK) XMJNK=XM1NK-1. ■- 40 

KXCYC*XMAXK-XMjNK 41 

IF (KXC YC.GT.7. AND.KCYMAX.EQ.7 ) GO TO 13 42 

IF (KXCYC.GT.15.AND.KCYMAX.EQ.15) GO TO 14 43 

SCALEX=SIZ£/FLOATIKXCYC) 44 

>XKIN=XKINK 45 

XXDEL*VAL(5)/SCALEX 46 

NCX=5*KXCYC 47 

NCXPTS=9*KXCYC 48 

CTH=THETA*RAD 49 

$TH*SIN (CTH I 50 

CTH=COS(CTH) 51 

C THY= ( THTY-90. ) *RAD 52 

STHY=SIN(CTHY) 53 

CTHY=COS(CTHY) 54 

DC 3 1=1,9 55 

XSPC(I)*TIC(I) *CTH* SCALEX 56 

2 YSPC<I)*TIC(1 1 *STH4 SCALEX 57 

XN=X+CTH*Sl ZE 58 
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YMY+STH*SIZE 

XNXaXN 

YNY=YN 

CALL PICT ( XNt YN*3 ) 

OC 4 1 = 1 »NOX PT$ 

JJ=I-1 

KK=tMOD( JJ t 9Hl 
XN=XN-XSPC(KK> 

YN*YN-YSPC(KK) 

A CALL SYMBOL LXN , YN » . 1 1 13 .THET A f -2 > 

DC 5 1*1 ,4 

XSPC< I ) = XSPC <2*I**1 l+XSPC(2M) 

5 yspc< n*YSPctz«r-i) ♦yspcc?*!! 

X$PC(5>*XSPCf9) 

Y$PC ( 5 ) * YSPC ( 9 ) 

XNU=XNX-STHY*HGT 

ym*yny+cthy*hgt 

KOTO*1 

IF (ABS(XSPC)*LT«1* E-3 ) KOTO*2 

OC 8 1=1 tNOX 

JJ*I-1 

KK=MOD( J J* 5 ) + l 
CALL PLOT <XNU,YNU,3) 

CALL PLOT (XNX,YNY,2> 

XN X = XNX-XSPC <KK ) 

YNY=YNY-YSPC ( KK ) 

XNl*XNU-XSPC<KK) 

YNU* YNU— YSPC ( KK > 

GC TO t 6 »7) « KOTO 

6 TEMM*YNU 
YM = YNY 
YNY=TEMM 
GC TO 8 

7 TEMM-XNU 
XNfXNX 
XNX-TEMP 

6 CONTINUE 

CALL PLOT ( 0* « O* »3 ) 

NCTE»N0X4:1 

DXB=-.05 

DYB=-*05*{n-3.*STHU$IGN<*15 t NNI 

DXC*S!7E/2«- f# 12*FLCAT H I ABS( NN)+7 ) /2 ) ) 

DYC*(-*075 + SIGNU375tNNM 

XN*X*-DXB*CTH-DYB*STH 

YN = Y+OYB*CTH*!DXB*$TH 

SPC*.20 

SPS*.06 

XN T* SORT ( DXC**2+DYC**2 ) 

DO 12 I *ltNCTE 
JJ»I-l 

KK=MOD< JJ*5> 

L L *5-KK 

IF (KK* NE# 0 ) GO TO 9 

CALL SYMBOL < XN tYN t «1 f 2H10«0« »2) 

X T *XN+SPC 
YT=YN+SPS 

CALL NUMBER CXTtYTf07tXMrNK»0.,-lI 
XM INK*XM INK+ 1* 

GC TO 10 

S CALL NUMBER I XN*YN # *1 f VAL ( LL > «0« t -ll 
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1C XN=XN+XSPC (LL) 

YN=YN+YSPC(LL> 

GC TO I 11,12) » KATE 

11 IF lXNT.GT.SQRT(XN**2*YN**2)) GC TO 12 
XT=X+OXC*CTH-DYC*STH 

YT = Y+DYC*CTH+D XC*STH 

CALL SYMBOL I XT» YT, .16 * BCD.THET A, I ABS t NN) ) 
KATE^2 

12 CONTINUE 
RETURN 

13 WRITE 16,17) 

GO TO 15 

16 WRITE 16,16) 

1* CONTINUE 

STOP 
C 

16 FORMAT (33HCNLY 15 CYCLES ALLOWED FOR X-AXIS) 

17 FORMAT (32HCNLY 7 CYCLES ALLOWED FOR Y-AXIS) 
END 
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Subroutine TRNSFM 


SUBROUTINE TRNSFM ( X,Y,N,THETX,THETY) 
DIMENSION XU), Y 1 1 ) 

RETURN 

END 
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Subroutine AXGRID 


C HIS SUBROUTINE WILL DRAW A LINEAR AXIS WITH ITS GRID LINES 1 

SUBROUTINE AXSRID I X, Y , BCD, N,SI ZE ,THET A , XMI N,DX,OV , HGT t THTY l 2 

3 

WF ERE- X,Y IS THE COORDINATE OF THE STARTING POINT OF THE AXIS. 4 

BOTH ARE FLOATING POINT AND PAGE INCHES. 5 

BCD IS THE LOCATION OF ALPHA INFORMATION FOR AXIS TITLE. 6 

USUALLY VARIABLE NAME. (NORMALLY SET UP WITH LITERAL) 7 

N IS THE NUMBER OF CHARACTERS IN BCD TITLE. A NEGATIVE 0 

N PLACES THE ANNOTATION ON THE CLOCKWISE SIDE OF AXIS 9 

LINE AND VICE-VERSA. 10 

SIZE IS THE LENGTH OF THE AXIS TO BE DRAWN. SIZE IS 11 

FLOATING POINT AND SHOULO BE MULTIPLY OF C 10.0/0 IV ) . 12 

THETA IS THE ANGLE OF THE AXIS MEASURED COUNTER-CLOCKWISE 13 

FROM X AXIS. THETA IS FLOATING POINT DEGREES. 14 

XMIN IS THE VALUE OF VARIABLE AT THE FIRST POINT OF THE 15 

AXIS. XMIN IS FLOATING POINT. (SEE NOTE) 16 

OX IS THE DIFFERENCE BETWEEN SECOND AND FIRST VALUE OF 17 

VARIABLE ALONG AXIS. DX IS FLOATING POINT. (SEE NOTE) 18 

DV IS THE NUMBER OF DIVISION PER INCH OF PAPER TO BE USED 19 

DV IS FLOATING POINT. (MAY BE 10.0,20.0.25.0,25.4) 20 

HGT IS THE LENGTH OF THE OTHER AXIS 21 

THTY IS THE ANGLE OF THE OTHER AXIS TO BE USED. 22 

NOTE- THE SECOND VERSION OF SCALE PLACES XMIN IN VAR I A8L E ( J*K + 1 ) AND 23 

DX IN VAR I ABLE ( J*K+K+1 ) • WHERE J IS NUMBER OF ELEMENTS IN 24 

ARRAY AND K IS THE REPEAT CYCLE OF MIXED ARRAY. 25 

26 

DIV*DV 27 

TWC=2.0 28 

IF (DIV.LT.1..0R.0IV.GT.25.4) DIV=10.0 29 

TH » THETA *0.0174532925 30 

THETY*(THTY-90. >*.0174532925 31 

STHYaSINITHETY) 32 

CTHY-COS 1 THETY ) 33 

C TH=C05 ( TH ) 34 

STH=SIN(TH) 35 

OIVA*OIV 36 

IF (0IV-20.C ) 1,2,2 37 

1 0 1 VA=t2.0*D I V 38 

Th0*1.0 39 

< SPACE=1C.0/C I VA 40 

DX8— .05*(1.+2.*CTH> 41 

0YB=-.O5*U.-5.*STHHSIGN(.15,N) 42 

DXC=SIZE/2.- (• 12*FLCAT ( (IABS(N) )/2) ) 43 

0YCM-.C75 + $IGN(.575,N) » 44 

DN0TE=SIZE/SPACE/2. 0+1.0 45 

NCTE^ONGTE 46 

IF I (ONCTE-FLQAT(NOTE) J.GT.0.99) N0TE=NCTE + 1 47 

SPC*CTH*SPACE*2.0 48 

SPS*STH*SPACE*2.0 49 

XN=X+DXB+CTH-DYB*STH 50 

YN = Y+DYB*CTF+’0X8*$TH 51 

ACY=ABS10X»*10.0/DIV 52 

EXsO.O 53 

IF (AOY) 3,7,3 54 

3 If (AOY-IOC.O) 6,4,4 55 

4 AOYaADY/lQ.C 56 

E> S E X+1.0 57 

GC TO 3 58 
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5 


ADY=ADY*10.C 59 

EX*EX-1.0 

6 IF ( ADY-1. 1 5,7*7 61 

7 AeSV=XMIN/lC.O«*EX 62 

NAO Y=AOY 63 

IF < ABS (AOY-FLOAT(NADY) ) »GT««99 ) ADY* 1 FIX ( ADY+S IGN( *01 » ADY ) ) 64 


ADY«SIGN<ADY,DX)*TWC 65 

N7*SQRT (DXC**2+DYC*»2) 66 

DC 10 1*1, NOTE 67 

CALL NUMBER ( XN »YN , .1 , ABSV,0. ,0 > 68 

AB SV*A6SV+A0 Y 69 

XN*XN+SPC 70 

YN*YfHSPS 71 

IF (NT) 10,8,10 72 

fi XT =X+OXC*CTH-DYC*STH 73 

YT*Y+DYC*CTH*OXC«<STH 74 

call Symbol i.xt, yt,.14,bcd, theta, iassind 75 

IF (EX) 9,10,9 76 

9 XT*XT+o 12*F10AT < I A85 t N) l*CTH 77 

YT*YT+ol2*FLOATUAB5(N) 1 *STH 78 

CALL SYMBOL ( XT,YT , .14 , 7H I X10 ), THETA, 7) 79 

XT=XT+(.12*5.*CTH-.06*$THI 80 

YT*YTM * 12*5s4STH+»06*CTH I 81 

CALL NUPBER (XT, YT,. 08, EX, THETA, -1) 82 

1C N7*NT-1 83 

DN7IC*DIV*SIZE/10. 0+1.0 84 

NT IC*DNT IC 85 

IF < C0NT1C-FLQAT(NTIC1).GT.0*99) NTIC=NTIC+1 86 

XN=X+CTH*SIZE 87 

YN*Y<-STH#SIZE 80 

XN X= XN 89 

YNY*YN 00 

XNU*XNX-STHY*HGT 91 

YNU=YNY*CTHY»HGT 02 

SPC=CTH*10.C/DIV 93 

SPS=STH *10*0/0 IV 99 

CALL PLOT (XN,YN,3) 95 

XN*FL0AT(NTIC-1)*SPC«-X 96 

YN*FLOAT«NTIC»X)*SPS+Y 97 

CO 11 1*1 ,NTIC 98 

CALL SYPBOL < XN ,YN,0.2*SPACE»13 ,THETA»-2 ) 99 

XN*XN~SPC 100 

11 YN*YN“SPS 101 

NTIC^NTIC-1 102 

KCTO^l 103 

IF 1 ABS (SPS5.LT.UE-03) K0T0-2 104 

DO 14 I =1 , NT I C 105 

CALL PLOT (XNU,YNUr3> 106 

CALL PLOT (XNU»YNU,3) 107 

CALL PLOT ( XNX »YNY , 2 1 108 

XNX=XNX— SPC 109 

VNY*YNY-SPS 11° 

XNli*XNU-SPC HI 

YMJ a YNU-SPS 112 

GO TO (12,13) , KOTO ll 3 

12 TEPH^XNli 114 

XNU*XNX 115 

XNX*TEPP 116 

GO TO 14 117 

13 T£PM=YNU 118 

YNL*YNY 1 * 

YNY*7EPP 1° 

14 CONTINUE 1^1 

RETURN I 22 

END I 23 " 
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APPENDIX E 


USER’S MANUAL 
HOW TO USE REDUCE H 

The input to REDUCE II consists of five basic pieces of information: 

(1) A card indicating whether the super G form is ARTIFICIAL or NATURAL 

(2) The block diagram inputs 

(3) The transfer function (or functions) desired 

(4) A set of algebraic expressions representing the block diagram 

(5) A sequence of variables separated by commas representing the reduction order 

string 

With the exception of the first card, the input data can be punched anywhere on a card. 


Examples of Input to REDUCE II 

Examples corresponding to the five basic types of input are described as follows: 

(1) This card must have ARTIFICIAL or NATURAL punched on it starting in card 
column 1. 

(2) The block diagram inputs are specified as, for example, 

X(7)$X(5) $X(1) $ $ 

A dollar sign separates each input. Two dollar signs in succession (no blanks between 
them) indicate the end of the block diagram inputs. There may be no more than five in- 
puts for a block diagram as the program is currently dimensioned, and only one card 
may be used. 

(3) The desired transfer function for the block diagram is given by, for example, 

X(2)$X(1) $ $ or X(2)$$ 

When two variables are indicated, REDUCE n solves the block diagram for the ratio of 
the first variable given (the output) to the second (an input). When one variable is given, 
REDUCE II solves for the ratio of this variable to all the inputs. Again two dollar signs 
in succession indicate the end of this information, and only one card may be used. 

(4) Block diagram expressions are written as 

(±N±G(A)*G(B)* . . . *G(F))*X(B)±G(M)#G(N)*. . . *X(C)±. . . $ 
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where A, B, C, ... are integers, N is a numerical constant, the G(i) s are functions, 
and the X(i)’s are variables. This allows very general forms that consist of variables 
X together with their coefficients, which are the sums and products of G’s and con- 
stants. Notice again that a dollar sign terminates an expression. No more than five 
cards may be used per expression. An example of an expression card is 

X(1)+X(2)-X(3)+G(4)*X(9)+(1. +G(1)*G(2))*X(4)$ 

FORTRAN notation is used to indicate subscripts and arithmetic operations. Also the 
=0 to make this a mathematical equation is implied. On the very last expression, the 
user must use two dollar signs in succession. This indicates to REDUCE II that there 
are no more expression cards to follow. 

(5) The reduction order string terminates the input to REDUCE II. Typically, it is 

X(A), X(B), X(C), , X(D), X(E), , X(F), . . . X(M)$ 

where the A, B, C, . . . are integers. The double commas partition the string into sub- 
strings and indicate a super G substitution. The dollar sign terminates the reduction 
string and should not be immediately preceded by any commas. An example is 

X(2), X(5), X(6), , X(10), X(7), , X(8) $ 

There must be no repetition of a variable in the reduction string. All variables but the 
output and inputs must be represented in the reduction string. Any number of cards may 
be used. Notice also that one dollar sign terminates the reduction string. These are all 
the data that REDUCE II needs. 


REDUCE II Output 

The output consists mainly of a set of cards. Included in these cards are the super 
G’s, the transfer function (or functions), and the transfer function identifiers. EVAL II 
uses this punched output in several places. The punched deck, after execution of 
REDUCE II, appears in the sequence shown in figure 17. 

REDUCE II prints a copy of the super G cards and the transfer function. Also the 
user's input in printed. Should REDUCE II detect an error, a message is printed and 
execution is usually stopped. A list of the messages and their causes follows: 

Message Cause 

(1) More than 600 variables The number of variables exceeds the maxi- 

mum allowed. 
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Message 

(2) More than 600 equations 

(3) More than 1200 functions 

(4) More than five block diagram inputs 

(5) More than five cards for a block 

diagram expression 

(6) The solution may be in error be- 

cause equation was not 

eliminated, equation is 

(7) Number of variables in order 

string is not consistent with 
the number of expressions - 
execution stopped 

(8) Only one equation containing X( ) 


(9) Warning - the number of super 
functions (super G's) exceeds 
the maximum number of functions 
(G T s) permitted, therefore sub- 
stitution process stopped 


(10) Subscript too large for this version 

(11) Invalid option specified on super G 

from card 


Cause 

The number of equations exceeds the maxi- 
mum allowed. 

The number of G functions before the reduc- 
tion starts exceeds the maximum allowed. 

The number of block diagram inputs exceeds 
the maximum allowed. 

Block diagram algebraic expressions may 
not use more than five cards. 

A variable used in the block diagram was not 
eliminated. This is because it was not in- 
cluded in the order string. 

There are probably one or more block dia- 
gram expressions missing, or there is a 
mistake in the order string. Check the 
input data. 

Only one equation was found containing the 
variable indicated. This indicates an 
error in the order string or in an equation. 
This error message will also be printed 
if there are no equations with the variable 
indicated. 

No more super G's will be created since the 
dimension of G is exceeded. Although the 
super G substitution process will stop, 
the reduction process will continue. It is 
likely that the expressions will become 
very large and that FORMAC storage 
space will be exceeded. 

Subroutine INSTRN has found a subscript 
larger than 600 in the reduction string. 

User has not specified NATURAL or 
ARTIFICIAL form. 
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4 


Transfer function 
identifier card 



HOW TO USE EVAL H 

Before EVAL II can be used, the output from REDUCE II must be rearranged. The 
FORTRAN-Iike cards go into subroutine COEFF of EVAL H. The other cards are data 
to identify transfer functions. They also tell EVAL II how many transfer functions are 
to be evaluated for any particular run. The cards shown in figure 17 are regrouped as 
shown in figure 18. 



These cards immediately 
follow the $DATA card. 


These cards all go into COEFF. 
The COMPLEX GU t DENOMO 
card goes ahead of all executable 
statements. The rest follow 
a comment card in COEFF. 

See listing of COEFF. 


Figure 18. - REDUCE II output regrouped. 
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The user’s functions or G’s should be placed before the super G cards. These 
functions may be of any form as long as they are in FORTRAN. Should they be one of 
the following forms, the namelist input section of EVAL II can be used: 


singly as 2 + bs + c SINHF(I) 

(El) 

cosh^y^as 2 + bs + c^ = COSHF(I) 

(E2) 

sin^dy^as 2 + bs + c^ = SINF(I) 

(E3) 

cos ^d y^as 2 + bs + c j = COSF(I) 

(E4) 

exp (d y^as 2 + bs + cj = EXPF(I) 

(E5) 

yas 2 + bs + c = SQRTF(I) 

(E6) 


where I is an integer constant. 


EVAL II Input 

Section I . - The first data cards to EVAL IT are the transfer function identifiers that 
REDUCE II punched out. 

Section II . - The next card of the EVAL input must contain $INPUT beginning in 
card column 2. This starts the namelist section of EVAL n data. The following input 
forms may appear in any order, but all must begin in card column 2: 

(1) NEXP = Number, EXPA = List, EXPB = List, EXPC = List, EXPD = List, 

(2) NSIN = Number, SINA = List, SINB = List, SINC = List, SIND = List, 

(3) NCOS = Number, COSA = List, COSB = List, COSC - List, COSD = List, 

(4) NSINH = Number, SINHA = List, SINHB = List, SINHC = List, SINHD = List, 

(5) NCOSH = Number, COSHA = List, COSHB = List, COSHC = List, COSHD = List, 

(6) NSQRT = Number, SQRTA = List, SQRTB = List, SQRTC = List, 

(7) NK = Number, K = List, 

(8) NT = Number, T = List, 

(9) NZ = Number, Z = List, 

(10) NWN = Number, WN = List, 
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(11) FSTART = Initial frequency, FEND = Final frequency, DELTAF = Increment 
(Frequencies are specified in hertz.) 

(12) LENGTH = Length of the abscissa in inches. 

Any of these forms may use more than one card as long as each card ends with a 
comma. The suffixes A, B, C, and D in the function names represent the variables 
implied in the function argument. (Note: no D inSQRT. ) For example, EXP A, EXPB, 
EXPC, and EXPD represent the variables a, b, c, and d in the function 


exp 



+ bs + c 


(See eqs. (El) to (E6) for additional forms. ) 

The term "Number, *» as used in the preceding forms, means an integer number 
which represents the number of a particular type of function or variable being used. For 
the functions, this number must not be greater than 5; for variables, not greater than 
10. The term "List" means a list of values that the variables are to take on. The ele- 
ments of this list must be separated by commas and must be in the proper order. For 
example, NK=2, K=l, 6. 25E-2 would set K(l) equal to 1.0 and K(2) equal to 0. 0625. 

For consecutive runs (runs where the transfer function form does not change), only the 
changed numbers need to be given on a new input. For example, 


$INPUT 
K(2)*l. 25E-2, 
$ 


A card with a $ in card column 2 must be the very last card of section n. This 
terminates the namelist data. 

Section III (optional) . - If LENGTH in the previous section were set to zero, no 
plots would be made. A positive LENGTH signals EVAL II to make CALCOMP plots. 
This means that EVAL H needs (1) a plot title and (2) information to make linear or log 
plots. This information the user punches onto two separate cards. 

(1) PLOT TITLE CARD - The user may use all 80 card columns for the plot title. 

(2) PLOT TYPE CARD - To get a set of semilogarithmic plots, the word LOGPLOT 
must appear on this card. LOGPLOT must start in card column 1. A blank card is suf- 
ficient to get linear plots. 

The user may run several cases of data in one execution of EVAL II. That is, many 
combinations of system information may be examined in sequence. To do this, the user 
may change any of the data values in section H. After completing the frequency band 
FSTART to FEND in steps of DELTAF, EVAL n returns to the statements that read the 
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section II (namelist) data. "EVAL II would - read- new data (or stop -execution if there is no 
data) and start anew. Remember that if LENGTH is left greater than zero, a new plot 
title and new plot type information must be supplied or the program will execute improp- 
erly. 


EVAL II Output 

The output printed by EVAL II lists the following items: (1) All the numerical input; 
(2) a table to show the correspondence between the transfer function and its identifying 
number; (3) the frequency band FSTART to FEND in DELTAF; and (4) the identifying 
number, frequency (Hz), real part, imaginary part, absolute value, phase angle (deg), 
and magnitude (dB) of the transfer function. 

The CALCOMP plots are of two types. The plots on linear axes have absolute value 
and phase angle as a function of frequency. Those on semilogarithmic axes have deci- 
bels as a function of frequency and phase angle as a function of frequency. Each plot is 
identified with the user-supplied title and the transfer function given as X^/Xg, where 
A and B are integers. 


A WORKED EXAMPLE 



The block diagram shown is quite simple but demonstrates the use of REDUCE II. 
We want to solve this diagram for Xj to X£ using naturally formed super G’s. The 
equations associated with the block diagram are 

Xj - G,x 2 = 0 
x 4 - GgXg = 0 
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X 1 - G 3 X 5 = 0 

X 2 - x i - x e = 0 


X 0 + X„ - X c = 0 


The input data to REDUCE II for this example would be 


NATURAL FORM 
X ( 2 ) $ $ 

xm$$ 

X ( 3) -G( 1 ) *X (2 ) * 

xm-xm-x<6>$ 

XUK(2)*XI6U 
XI 3 ) +X{ A) -X ( 5 ) $ 
X(1)-G(3>*X(5)$$ 

X( 3) » X( A) t * X( 6 ) , X ( 5) $ 


The printed output from REDUCE n for this example is 

THE BLOCK D IA GRAM INPUTS 
X(2H 

the SCUJTION IS FOR Xlllt TO THE INPUTS 

the block diagram expressions 

1 X ( 3 )—C- ( l )*X ( 2 It 

2 X(21-X( Il-XfCJ* 

3 X { 4 ) — G ( 2 I*X ( 6 )S 
A X ( 3 ) +X ( A )— X ( 5 H 

5 xu >-cm*xi ' 

THE REDUCTION ORDER IS 

l X(3) v X(4)tf X(6)fXm$ 

THE LARGEST VARIABLE SUBSCRIPT IS 6 

THE NUMBER OF EXPRESSIONS IS 5 

THE LARGEST FLNCTION (0) SUBSCRIPT IS 3 


* * * OUTPUT FROM REDUCE II * * * 

TRANSFER FUNCTION X( 1I/X< 21 

THE NUMERATOR 

G! 1)*GC3)+G(2)*G(3> 


THE DENOMINATOR 
Gf 2)*GC3) + 1.0 



The cards punched out by REDUCE II are 


DC* 

*G { 2 ) *G( 3 ) +1 ■ 0 
NC( H = 

*Gm*Gm+G<2) *G<3> 

CCNPLEX G< 3 f. >DENOM ( i) 

1 

1 

2 


The FORTRAN expressions for Gj, Gg, and Gg are 

G(1)=-K(1)*Z(1)*Z(1)*Z(2)*S*S*SINHF(1)/ 

(T(I)*Z(3)*Z(4)*SQRTF(l)*COSHF(l)) 

G(2)=K(1)+Z(5)*S 

G(3)=1./(Z(6)*S*S) 

These expressions immediately precede the super G output in COEFF. See page 72 for 
a typical example. The constants associated with this example are 

Kj = 7.6336x10 s 

Tj - 1.0416X10' 2 


Zj = 0.75 
Z 2 = 62.4 
Z 3 = 6.48X10 7 
Z 4 - 32.2 
z 5 = 10 
z 6 = i 


SINH(1) 

COSH(l) 


-SINK 


|d • "^c + bs + j 
COS ll/d • + bs + 


as 
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where 


a = 4. 1006x10 

b = 10“ 4 
c — 0 
d = 3.333 

The following cards are the input to EVAL II for this example: 


1 

1 

2 

$ inpit 

NSINMlt SINHA=4*IOO625E-6.SINHB*l*OE-4*SINHC»O..SlNH0»3*3333 
NC0SM1* C0SHA»4* 1006256-6 » C0SHB*1 *0E-4 * COSHC*0* »C0SH0 st 3*3333 
N SQR 1 *1 . SORT A*4*l00625E-6* SQRTB*l* E-4 , SQRTC-0. * 

NK-1 *K^«76335936E6 * 

N>T*1»T* 1*04166-2* 

NZ=6 * Z** 750*62*4 »6*48E 7 *32*2*10* *1* * 

FSTABT*l*t*FENO»50* *DELTAF**1* * LENGTH 3 10 * * 

& 

A WOftBeC EXAMPLE 
LINEAR PLOT 


The output from EVAL II for this example is 



Fortran evaluation of formac output 

CALCOWR PLOTS HAVE BEEN REQUESTED 
A WORKED EXAMPLE 


SINHA 

0. 41006E-05 
COSHA 

O.A130&E-05 
SORT A 

0.41D06E-05 

K 

0.76336E+D6 


T 

0. 104166-01 


I l 

1 0. 75 OOO E+OO 

2 0.62400E+02 

3 0.64800E+08 

4- 0.32200E+Q2 

5 0* 1330DE+02 

6 O.lOOOOE+ai 


S1NH8 S I VHC SINHD 

0.10000E-03 0. 0.333336*01 

CDSHB COSHC CDSHD 

0. 10000 E-03 0. 0.333336*01 

SQRTB SORTC 

0*10000 E-03 0* 


NUMBER TRANSFER FUNCTION 

1 X( 1J/X4 21 


FSTART* 

NUMBER 


1.00000 

FEND* 

50.00000 

DEL TAF= 1 

.00000 

REQUENCV 


REAL PART 

IMAGINARY PART 

ABSOLUTE VALUE 

1.000000 


1.000264 

-0.000001 

1.000264 

2.000000 


1.001059 

-0.000004 

1.001059 

3.000000 


1.002390 

-0.000014 

1,002390 

4.000000 


1.004264 

-0.000034 

1.004264 

5.000000 


1.006695 

-0.000067 

1.006695 

6.000000 


1.009698 

-0.000118 

1.009698 

7.000000 


1.013295 

-0.000190 

1.013295 

8.000000 


1.017510 

-0.000291 

1,017510 

9.000000 


1.022374 

-0.000425 

1.022374 

lo.oonooo 


t ,027924 

-0.000600 

1.027924 

11.000000 


1.034204 

-0.000825 

1.034204 

12.D0DD00 


1.041265 

-0.001112 

i. 041266 

13.000000 


1.049170 

-0.001472 

1.049172 

14.000000 


1.057993 

-0.001923 

1.057995 

15.000000 


1.067821 

-0.002485 

1.067824 

16.000000 


1.078758 

-0.003105 

1.078763 

17.000000 


1.090931 

-0.0040 54 

1.090936 

18.000000 


1.104491 

-0.005135 

1.104503 

19.000000 


1.119624 

-0.006485 

1.119643 

20.000000 


1.136558 

-0.008175 

l. 136587 

21.000000 


1.155573 

-0,010304 

1.155619 

22.000000 


1. 177022 

-0.013006 

1.177094 


PHA5E ANGLE 
-0,000030 
-0.000238 
-0.000807 
-0.001930 
-0.003810 
-0.006671 
-0.010762 
-0.01636* 
-0.023798 
-0.033*37 
-0.045719 
-0.061162 
-O.0BO388 
-0.104148 
-0.133356 
-0.169141 
-0.212902 
-0,266397 
-0 ,331851 
-0.412119 
-0.510904 
-0.633071 


AMPUTIDE IN 08. 
0.002296 
0.009195 
0.020730 
€.036958 
0.057957 
0.003832 
0.114717 
0.150772 
0.192196 
0.239220 
0.292123 
0.351 231 
0.416930 

0.489471 
0.565991 
0.658519 
0,756004 
0.863339 
0.9B1593 
1.1 12056 
1.256293 
1.416223 


1 23.000000 

l 24.000000 

1 25.000000 

1 26.000000 

1 27.000000 

1 28.000000 

l 29.000000 

1 30.000000 

1 31.000000 

1 32.000000 

1 33.000000 

1 34,000000 

1 35,000000 

1 36.000000 

1 37.000000 

l 38.000000 

1 39.000000 

1 40.000000 

1 41.000000 

1 42.000000 

1 43.000000 

1 44.000000 

1 45.000000 

1 46.000000 

1 47.000000 

l 48,000000 

1 49.000000 

1 50.000000 

PLOT MADE 
PLOT MADE 


1.201352 

1.229135 

1.261115 

1.29B270 

1.341946 

1,393926 

1.456726 

1.533872 

1.630300 

1.752605 

1.907810 

2.094135 

2.254446 

2.110127 

1.178423 

0.153711 

-0.136001 

-0.071503 

0.061374 

0.166983 

0.292803 

0.3B0371 

0.453432 

0.515370 

0.568820 

0,615761 

0.657672 

0.695668 


-0.016463 

-0.020935 

-0.026797 

-0.034801 

-0.045184 

-0.059864 

-0.080787 

-0.111610 

-0.158890 

-0.235092 
-0.365391 
—0.602963 
-1.053980 
-1.819422 
-2.424559 
-2.014335 
-1.303576 
-0.836700 
-0.569645 
-0.412H9 
-0.313532 
-0.2 48257 
-0.202922 
-0.170171 
-0.145730 
-0.126993 
-0.112305 
-0.100572 


1.201465 

1.229313 

1.261400 

1.298739 

1.342707 

1.395211 

1,458964 

1.537927 

1.636024 

1.76B303 

1.942465 

2.179212 

2.480654 

2.786204 

2.695766 

2.020292 

1.310651 

0.839758 

0.572942 

0.452553 

0.428994 

0.454217 

0.496768 

0.542738 

0.567191 

0.626721 

0.667192 

0.702900 


-0.785106 
-0.975B05 
-1.217290 
-1.526637 
-1.92B449 
-2.459133 
-3.174252 
-4.161703 
-5.566511 
-7.639974 
-10.842202 
-16.062666 
-25.056649 
-40.769007 
-64.078579 
-65,636287 
-95.956075 
-94.884479 
-83.650584 
-65.595722 
-46,958082 
-33.1312T2 
-24.1097H 
-1 B . 27282 4 
-14.369910 
-11.653161 
-9.69D418 
-8.226227 


1.59422Q 
1.793248 
2.017056 
2.2 70*40 
2.559623 
2.89 2 795 
3.280894 
2.738714 
4,286406 
4.951 131 
5.767153 
6.765591 
7.919291 
6. 900259 
8.6136S1 
6.10 7851 
2.349741 
-1.516920 
-4.837784 
-6.066604 
-7.35097® 
-6.054729 
-6.076932 
-5.308194 
-4.624406 
-4.030047 
-3.514981 
-3.062 131 
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